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ABSTRACT 

Herbig Ae/Be stars span a key mass range that links low and high mass stars, 
and thus provide an ideal window from which to explore their formation. This paper 
presents VLT/X-Shooter spectra of 91 Herbig Ae/Be stars, HAeBes; the largest spec¬ 
troscopic study of HAeBe accretion to date. A homogeneous approach to determining 
stellar parameters is undertaken for the majority of the sample. Measurements of the 
ultra-violet (UV) are modelled within the context of magnetospheric accretion, allow¬ 
ing a direct determination of mass accretion rates. Multiple correlations are observed 
across the sample between accretion and stellar properties: the youngest and often 
most massive stars are the strongest accretors, and there is an almost 1:1 relationship 
between the accretion luminosity and stellar luminosity. Despite these overall trends 
of increased accretion rates in HAeBes when compared to classical T Tauri stars, we 
also find noticeable differences in correlations when considering the Herbig Ae and 
Herbig Be subsets. This, combined with the difficulty in applying a magnetospheric 
accretion model to some of the Herbig Be stars, could suggest that another form of 
accretion may be occurring within the Herbig Be mass range. 

Key words: stars: early-type - stars: variables: Herbig Ae/Be - stars: pre-main 
sequence - stars: formation - accretion, accretion discs - techniques: spectroscopic. 


1 INTRODUCTION 

Herbig Ae/Be stars are pre-main sequence, PMS, stars that 
bridge the key mass range of 2-10 Mq , between low and high 
mass stars. Their importance lies in linking the reasonably 
well understood formation of low mass stars, like the numer¬ 
ous Classical T Tauri stars, CTTs (which have M* < 2 M© 
and will go on to form stars like our own Sun), to the rarer, 
more deeply embedded high mass stars (or massive young 
stellar objects, MYSOs). The greater number of HAeBes 
compared to forming high mass stars, along with them be¬ 
ing closer and optically visible, makes them a powerful link 
and important tool in furthering our understanding of star 
formation upto the high mass star format ion regime. 

HAeBes were originally identified by iHerbid (Il96d ') in 
an attempt to push the mass boundaries in understand¬ 
ing of CTTs, and had to meet the criteria of: “spectral 
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type A or earlier with emission lines; lies in an obscured 
region; the star illuminates fairly bright nebulosity in its 
immediate vicinity”. The latter two criteria have been re¬ 
laxed since then in order to find more potential targets 
(see iFinkenzeller fc Mundd[l984l : iThe. de Winter fc Perej 
1 1994 Vieira et aLlkOOsh . In these surveys more attention 
has been drawn towards the colours of the objects, par¬ 
ticularly in the infra-red, IR. This is because the IR often 
displays an excess of emission, compared to main-sequence, 
MS, stars of the same spectral type, and is the product 
of the circumstellar disc arou nd HAeBes. This has been 
confirmed in nume rous studies l|v3-^i den Ancker et al. I bOOd: 
iMeeus et al. 200l|), and by direct observations in th e op- 
tical llMcCaughrean fc 0’delllll99d : lGradv et al.lf200lh . sub- 
mm llMmmhigsfc Sajgenl| j^ and of scattered, polarised 
light ( Vink et al. 2002 , 2005l l . It is these indicators which 
ultimately point to the stars being young and in the PMS 
phase. T his young nature, s uspected by Herbig, was con¬ 
firmed bv IStrom et al.l (Il972l l who observed the HAeBes to 
have lower surface gravities than MS stars. With their PMS 
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nature established, the next obvious questions are: what 
happens with the disc-star interaction; how do they both 
evolve; and what physical mechanisms are present in these 
interactions? And of course, how are these aspects related 
to each other between the CTTs and massive young stellar 
objects? 


With regards to the disc-star interaction, an 
ultra-violet, UV, ex cess in the CTTs ijOarrisonl Il978l : 
ICullbring et al.lll998h has been shown to match well with 
the theory of disc-to-star accre tion within the Magneto- 
spheric Accretion, MA , regime (ICalvet fc Gullbrind 119981 : 
ICnllbring et al.l l200fll : llnglebv et al. 20 bil l. Under this 
paradigm the disc is truncated by the stellar magnetic 
held lines, and from here the material is funnelled by 
the held lines, in free-fall, onto the star. The accreted 
material shocks the photosphere causing X-ray emission, 
the majority of which is then absorbed by the surroundings, 
heating them, and is re-emitted at longer wavelengths; 
givin g rise to an observable UV excess (ICalvet fc Gullbringl 
Il998h . Therefore, measurement of the UV excess can be 
directly related to accretion from the disc to the star. The 
accretion rate has large implications on PMS systems; it 
affects the achievable hnal mass of not just the star, but 
also any possible planets forming with in the disc and the 
timescales upon whi ch they can evolve dLubow 8z D’Angelo! 
120061 : IPunhilll l201fih . Establishing an accretion rate also 
requires accurate stellar parameters which, until now, have 
mostly been performed on an ad-hoc basis dependent on 
the particular study. The results of this work will also 
help future studies to disentangle the complexities in the 
environments of these such as: outflows, infalling material, 
structure, and even possible ongoing planet formation. 

However, for MA to be applicable the star must have a 
magnetic field with sufficient strength to tr uncate the disc; 
for CTTs this i s of the order of kilo-Ga uss llGfioslLfcLamd 


19791: 


Boiivier et al 


Koenig] Il99ll : IShu et akl 1 19941 : Ijohns-Krulll 1200^ 
l.| 2007^ . For MS stars, magnetic fields driven 


by convection are not predicte d to exist for stars with 
Teff > 8300 K dSimon et al.ll2002ll . However, this may not be 
the case for PMS stars as th ere have been a few detections of 
magnetic fields in HAeB es l|Wade et al.ll200^ : ICatala et al.l 
120071 : iHubrig et al.ll200^ . but the origin of their fields re¬ 
mains unknown (be they dynamo generated or the result 
of fossil fields). The largest survey in to the magnet i c field s 
of HAeBes was performed recently bv lAlecian et al.l (l2013ll . 
and yields clear detections of only 5 stars out of 70. Wh en 
applying the theory of MA (lKoenigll99ll : [^hu et al.ll994l l to 
HAeBes, a weak dipole magnetic field of only a few hundred 
gauss, or even less, is needed for MA to occur (jWade et al.l 
I2OO7I : Icaulev fc Johns-Krullll2014| j. These strengths are be¬ 
low current detection limits, meaning that MA acting in 
HAeBes is still a possibility. 

A key aim of this paper is to provide the largest 
survey on direct accretion tracers in HAeBes to date. 
To do this, measurements of the UV excess are made 
and fitted within the context of MA shock modelling. 
This method of accretion-shock modelling has been suc¬ 
cessfully adapted and applied, in the m ajority of cases, 


2 OO 4 I: Donehew & Br 

ittainl 2 OIII: Mendiffutia et alJ l2011bl: 

Pogodin et al.l 20121: 

Mendigutfa et al. l2013l. 2 OI 4 I 1 . How- 


ever, the number of HBes analysed in previous works is of¬ 


ten small, particularly for early-type HBes, and needs to be 
tested further. The derivation of accretion rates, and other 
properties, depend heavily upon basic stellar parameters. 
Thus far, most works on accretion in HAeBes build upon 
previous stellar parameter determinations from a wide range 
of sources and methodologies. This can result in parameters 
which are not directly comparable within a sample. There¬ 
fore, our approach is to provide a homogeneous determi¬ 
nation of stellar parameters for as many of the targets as 
possible. This not only helps in our accretion determination 
but also helps in future works which require basic stellar 
parameters i.e. detailed modelling of the circumstellar disc, 
or energy balance in the SED. 

The overall aim of this paper is to provide a quantita¬ 
tive look into the properties of HAeBes, with a particular 
emphasis on how the accretion rate varies as a function of 
stellar parameters; along with an assessment of the applica¬ 
bility of using MA to obtain the accretion rate. To do this 
we present 91 HAeBe objects observed with the X-Shooter 
spectrograph, VLT, Chile. The paper is broken down into 
the following sections: Section [2] details the sample, obser¬ 
vations, and data reduction; Section [3] details the methods, 
and results, of deriving the basic stellar parameters; Sec¬ 
tion |4] presents the methods of measuring the UV excess; 
Section [S] presents the derived accretion rates, including a 
detailed description of how MA is applied to the HAeBe 
sample; Section |6] forms the discussion; and Section [3 pro¬ 
vides the conclusions of this work. Finally, photometric data 
and literature information on the sample are provided in the 
appendix. 


2 OBSERVATIONS AND DATA REDUCTION 


2.1 X-Shooter and Target Selection 

Observations were performed over a period of 6 months be¬ 
tween October 2009 and April 2010 using the X-Shooter 
echelle spectrograph - mo unted at the VLT, Cerro Paranal, 
Chile llVernet et al.ll201lh . X-Shooter provides spectra cov¬ 
ering a large wavelength range of 3000-23000 A, split into 
three arms and taken simultaneously. The arms are split 
into the following: the UVB arm, 3000-5600 A; the VIS arm, 
5500-10200 A; and the NIR arm, 10200-24800 A. The small¬ 
est slit widths available of 0.5”, 0.4” and 0.4” were used to 
provide the highest possible resolutions of R ~ 10000, 18000 
and 10500 for the respective UVB, VIS and NIR arms. In to¬ 
tal 91 science targets were observed in nodding mode using 
an ABBA sequence. Table [1] includes details of each targets 
RA and DEC, exposure times, and the signal-to-noise in 
each arm. The SNR is calculated by analysing a 30 A region 
of spectra centred about the wavelengths of 4600, 6750, and 
16265 A for the UVB, VIS, and NIR arm respectively. These 
regions were chosen as they are generally the flattest con¬ 
tinuum regions in each star. Although emission lines, and 
absorption lines of cooler objects, can artificially lower the 
measured SNR, for a fair treatment we stick with the above 
regions to provide a rough guide to the quality of each spec¬ 
trum. 

The targets were select ed fr om the catalogues of 
iThe. de Winter fc Pered I^1994^. and I Vieira et al.l ll 2003ll. 5 1 
targets were selected from The, de Winter fc Perez lll994h . 
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Table 1. Column 1 shows the target names, columns 2 and 3 are RA and DEC, column 4 gives the observation date, columns 5-7 give the exposure 
times for each arm, column 8 is the number of Detector Integration Times for the NIR arm, and finally columns 9-11 give the Signal-to-Noise 
Ratio in each arm. 


Name 

RA DEC 

(J2000) 

Obs Date 
(yyyy/mm/dd) 


UVB 

Exposure Time (s) 
VIS 

NIR 

NDIT 

UVB 

SNR 

VIS 

NIR 

UX Ori 

05:04:29.9 

-03:47:16.8 

2009-10-05 


90 

X 

4 


90 

X 

4 


5 

X 

4 

18 

42 

257 

306 

PDS 174 

05:06:55.4 

-03:21:16.0 

2009-10-05 


300 

X 

4 


300 

X 

4 

(20 

X 6) 

X 

4 

2 

127 

64 

225 

V1012 Ori 

05:11:36.5 

-02:22:51.1 

2009-10-05 


250 

X 

4 

(125 

X 2) 

X 

4 


20 

X 

4 

12 

128 

157 

229 

HD 34282 

05:16:00.4 

-09:48:38.5 

2009-12-06 

(15 

X 2) 

X 

4 

(15 

X 2) 

X 

4 


10 

X 

4 

3 

171 

215 

170 

HD 287823 

05:24:08.1 

02:27:44.4 

2009-12-06 

(15 

X 2) 

X 

4 

(15 

X 2) 

X 

4 


15 

X 

4 

2 

99 

171 

256 

HD 287841 

05:24:42.8 

01:43:45.4 

2009-12-06 


90 

X 

4 


90 

X 

4 


15 

X 

4 

6 

78 

200 

323 

HD 290409 

05:27:05.3 

00:25:04.9 

2010-01-02 

(15 

X 2) 

X 

4 

(15 

X 2) 

X 

4 


15 

X 

4 

2 

195 

157 

130 

HD 35929 

05:27:42.6 

-08:19:40.8 

2009-12-17 

(10 

X 2) 

X 

4 

(10 

X 2) 

X 

4 


5 

X 

4 

8 

36 

125 

161 

HD 290500 

05:29:48.0 

-00:23:45.8 

2009-12-17 


150 

X 

4 


150 

X 

4 


75 

X 

4 

2 

208 

341 

331 

HD 244314 

05:30:18.9 

11:20:18.2 

2010-01-02 

(15 

X 2) 

X 

4 

(15 

X 2) 

X 

4 


15 

X 

4 

2 

64 

140 

180 

HK Ori 

05:31:28.1 

12:09:07.6 

2009-12-17 


250 

X 

4 


250 

X 

4 


20 

X 

4 

12 

49 

120 

192 

HD 244604 

05:31:57.3 

11:17:38.8 

2009-12-17 

(15 

X 2) 

X 

4 

(15 

X 2) 

X 

4 


5 

X 

4 

8 

76 

170 

213 

UY Ori 

05:32:00.4 

-04:55:54.6 

2009-12-26 


300 

X 

4 


300 

X 

4 

(50 

X 6) 

X 

4 

1 

229 

326 

143 

HD 245185 

05:35:09.7 

10:01:49.9 

2009-12-17 

(15 

X 2) 

X 

4 

(15 

X 2) 

X 

4 


15 

X 

4 

2 

240 

171 

213 

T Ori 

05:35:50.6 

-05:28:36.9 

2009-12-17 

(25 

X 2) 

X 

4 

(25 

X 2) 

X 

4 

(2 

X 5) 

X 

4 

5 

156 

198 

132 

V380 Ori 

05:36:25.5 

-06:42:58.9 

2009-12-17 

(80 

X 2) 

X 

4 

(80 

X 2) 

X 

4 

(3 

X 2) 

X 

4 

20 

9 

178 

161 

HD 37258 

05:36:59.1 

-06:09:17.9 

2010-01-02 

(15 

X 2) 

X 

4 

(15 

X 2) 

X 

4 


10 

X 

4 

3 

62 

219 

168 

HD 290770 

05:37:02.5 

-01:37:21.3 

2009-12-26 

(15 

X 2) 

X 

4 

(15 

X 2) 

X 

4 


7 

X 

4 

4 

261 

244 

262 

BF Ori 

05:37:13.2 

-06:35:03.3 

2010-01-02 


90 

X 

4 


90 

X 

4 


10 

X 

4 

9 

25 

185 

329 

HD 37357 

05:37:47.2 

-06:42:31.7 

2010-02-05 

(10 

X 2) 

X 

4 

(10 

X 2) 

X 

4 


10 

X 

4 

4 

123 

209 

161 

HD 290764 

05:38:05.3 

-01:15:22.2 

2009-12-26 

(30 

X 2) 

X 

4 

(30 

X 2) 

X 

4 


7 

X 

4 

8 

61 

199 

233 

HD 37411 

05:38:14.6 

-05:25:14.4 

2010-02-05 

(25 

X 2) 

X 

4 

(25 

X 2) 

X 

4 


15 

X 

4 

3 

215 

183 

138 

V599 Ori 

05:38:58.4 

-07:16:49.2 

2010-01-06 


360 

X 

4 

(180 

X 2) 

X 

4 

(10 

X 2) 

X 

4 

10 

61 

149 

328 

V350 Ori 

05:40:11.9 

-09:42:12.2 

2010-02-05 


150 

X 

4 


150 

X 

4 


25 

X 

4 

6 

77 

260 

150 

HD 250550 

06:01:59.9 

16:30:53.4 

2010-01-02 

(15 

X 2) 

X 

4 

(15 

X 2) 

X 

4 


3 

X 

4 

10 

103 

218 

200 

V791 Mon 

06:02:15.0 

-10:01:01.4 

2010-02-24 


90 

X 

4 


90 

X 

4 


15 

X 

4 

6 

67 

343 

147 

PDS 124 

06:06:58.5 

-05:55:09.2 

2010-02-10 


300 

X 

4 


300 

X 

4 

(50 

X 6) 

X 

4 

1 

188 

62 

119 

LkHa 339 

06:10:57.7 

-06:14:41.8 

2010-01-17 


300 

X 

4 


300 

X 

4 


60 

X 

4 

5 

122 

63 

297 

VY Mon 

06:31:06.8 

10:26:02.9 

2010-02-08 

(150 

X 2) 

X 

4 

(100 

X 3) 

X 

4 

(2 

X 6) 

X 

4 

20 

11 

49 

102 

R Mon 

06:39:10.0 

08:44:08.2 

2010-02-01 


300 

X 

4 

(150 

X 2) 

X 

4 

(2 

X 6) 

X 

4 

20 

11 

49 

141 

V590 Mon 

06:40:44.7 

09:47:59.7 

2010-02-01 


300 

X 

4 


300 

X 

4 


50 

X 

4 

6 

143 

166 

288 

PDS 24 

06:48:41.8 

-16:48:06.0 

2009-12-16 


300 

X 

4 


300 

X 

4 


90 

X 

4 

3 

145 

48 

260 

PDS 130 

06:49:58.7 

-07:38:52.1 

2009-12-16 


300 

X 

4 


300 

X 

4 


60 

X 

4 

5 

125 

122 

273 

PDS 229N 

06:55:40.1 

-03:09:53.1 

2010-02-10 


300 

X 

4 


300 

X 

4 


100 

X 

4 

3 

111 

105 

191 

GU CMa 

07:01:49.6 

-11:18:03.9 

2009-12-16 

(2 

X 3) 

X 

4 

(2 

X 3) 

X 

4 

(2 

X 2) 

X 

4 

6 

127 

180 

142 

HT CMa 

07:02:42.7 

-11:26:12.3 

2010-01-30 


300 

X 

4 


300 

X 

4 


30 

X 

4 

10 

201 

245 

274 

Z CMa 

07:03:43.2 

-11:33:06.7 

2010-02-24 

(75 

X 2) 

X 

4 

(10 

X 3) 

X 

4 

(0.665 

X 5) 

X 

4 

20 

28 

87 

102 

HU CMa 

07:04:06.8 

-11:26:08.0 

2010-01-17 


300 

X 

4 


300 

X 

4 


50 

X 

4 

6 

185 

218 

217 

HD 53367 

07:04:25.6 

-10:27:15.8 

2010-02-24 

(2 

X 3) 

X 

4 

(2 

X 3) 

X 

4 

(2 

X 2) 

X 

4 

6 

62 

166 

220 

PDS 241 

07:08:38.8 

-04:19:07.0 

2009-12-21 


300 

X 

4 


300 

X 

4 


100 

X 

4 

3 

228 

295 

250 

NX Pup 

07:19:28.4 

-44:35:08.8 

2010-02-01 


120 

X 

4 

(60 

X 2) 

X 

4 

(2 

X 2) 

X 

4 

20 

44 

191 

125 

PDS 27 

07:19:36.1 

-17:39:17.9 

2010-02-24 


300 

X 

4 


300 

X 

4 

(2 

X 6) 

X 

4 

20 

9 

118 

228 

PDS 133 

07:25:05.1 

-25:45:49.1 

2010-02-24 


300 

X 

4 


300 

X 

4 

(40 

X 4) 

X 

4 

2 

3 

22 

50 

HD 59319 

07:28:36.9 

-21:57:48.4 

2010-02-24 

(10 

X 2) 

X 

4 

(10 

X 2) 

X 

4 


10 

X 

4 

4 

265 

14 

245 

PDS 134 

07:32:26.8 

-21:55:35.3 

2010-02-24 


300 

X 

4 


300 

X 

4 


150 

X 

4 

2 

196 

210 

157 

HD 68695 

08:11:44.3 

-44:05:07.5 

2009-12-21 

(20 

X 2) 

X 

4 

(20 

X 2) 

X 

4 


15 

X 

4 

3 

159 

179 

314 

HD 72106 

08:29:35.0 

-38:36:18.5 

2009-12-19 

(10 

X 2) 

X 

4 

(10 

X 2) 

X 

4 


10 

X 

4 

4 

49 

136 

154 

TYC 8581-2002-1 

08:44:23.5 

-59:56:55.8 

2009-12-21 


150 

X 

4 


150 

X 

4 


50 

X 

4 

3 

181 

205 

224 

PDS 33 

08:48:45.4 

-40:48:20.1 

2009-12-21 


300 

X 

4 


300 

X 

4 


150 

X 

4 

2 

278 

181 

186 

HD 76534 

08:55:08.8 

-43:27:57.3 

2010-01-30 

(10 

X 3) 

X 

4 

(15 

X 3) 

X 

4 


15 

X 

4 

3 

149 

259 

150 

PDS 281 

08:55:45.9 

-44:25:11.4 

2009-12-21 

(15 

X 2) 

X 

4 

(15 

X 2) 

X 

4 


7.5 

X 

4 

4 

193 

174 

203 

PDS 286 

09:05:59.9 

-47:18:55.2 

2009-12-21 


300 

X 

4 

(150 

X 2) 

X 

4 

(2 

X 6) 

X 

4 

20 

96 

205 

157 

PDS 297 

09:42:40.0 

-56:15:32.2 

2010-01-04 


300 

X 

4 


300 

X 

4 

(150 

X 2) 

X 

4 

2 

213 

210 

188 

HD 85567 

09:50:28.3 

-60:57:59.5 

2010-03-06 

(10 

X 3) 

X 

4 

(15 

X 3) 

X 

4 


2 

X 

4 

20 

142 

183 

283 

HD 87403 

10:02:51.3 

-59:16:52.7 

2010-03-06 

(15 

X 2) 

X 

4 

(15 

X 2) 

X 

4 


15 

X 

4 

2 

137 

137 

160 

PDS 37 

10:10:00.3 

-57:02:04.4 

2010-03-31 


300 

X 

4 


300 

X 

4 

(3 

X 6) 

X 

4 

15 

18 

182 

166 

HD 305298 

10:33:05.0 

-60:19:48.6 

2010-03-31 


90 

X 

4 


90 

X 

4 


45 

X 

4 

2 

155 

209 

207 

HD 94509 

10:53:27.2 

-58:25:21.4 

2010-02-05 

(20 

X 2) 

X 

4 

(20 

X 2) 

X 

4 


15 

X 

4 

3 

12 

171 

208 

HD 95881 

11:01:57.1 

-71:30:46.9 

2010-01-04 

(10 

X 3) 

X 

4 

(15 

X 3) 

X 

4 


2 

X 

4 

20 

56 

223 

296 

HD 96042 

11:03:40.6 

-59:25:55.9 

2010-02-05 

(10 

X 2) 

X 

4 

(10 

X 2) 

X 

4 


10 

X 

4 

4 

97 

107 

176 
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Table 1 — continued 


Name 

RA DEC 

(J2000) 

Obs Date 

(yyyy/mm/dd) 


UVB 

Exposure Time (s) 
VIS 

NIR 

NDIT 

UVB 

SNR 

VIS 

NIR 

HD 97048 

11:08:03.0 

-77:39:16.0 

2010-02-05 

(15 

X 3) 

X 

4 

(15 

X 3) X 

4 


2 

X 

4 

20 

208 

254 

200 

HD 98922 

11:22:31.5 

-53:22:09.0 

2010-03-30 

(3 

X 3) 

X 

4 

(3 

X 3) X 

4 

(0.75 

X 2) 

X 

4 

20 

85 

196 

103 

HD 100453 

11:33:05.3 

-54:19:26.1 

2010-03-29 

(10 

X 3) 

X 

4 

(15 

X 3) X 

4 


2 

X 

4 

20 

35 

110 

289 

HD 100546 

11:33:25.1 

-70:11:39.6 

2010-03-30 

(3 

X 3) 

X 

4 

(3 

X 3) X 

4 

(1 

X 2) 

X 

4 

20 

283 

206 

348 

HD 101412 

11:39:44.3 

-60:10:25.1 

2010-03-30 

(15 

X 2) 

X 

4 

(15 

X 2) X 

4 


10 

X 

4 

3 

89 

125 

169 

PDS 344 

11:40:32.8 

-64:32:03.0 

2010-03-31 


300 

X 

4 


300 X 

4 


150 

X 

4 

2 

236 

190 

157 

HD 104237 

12:00:04.8 

-78:11:31.9 

2010-03-30 

(3 

X 3) 

X 

4 

(3 

X 3) X 

4 

(0.75 

X 2) 

X 

4 

20 

25 

88 

216 

V1028 Cen 

13:01:17.6 

-48:53:17.0 

2010-03-29 


90 

X 

4 


90 X 

4 


10 

X 

4 

9 

101 

165 

148 

PDS 361S 

13:03:21.6 

-62:13:23.5 

2010-03-31 


300 

X 

4 


300 X 

4 


150 

X 

4 

2 

118 

192 

181 

HD 114981 

13:14:40.4 

-38:39:05.0 

2010-03-29 

(3 

X 3) 

X 

4 

(3 

X 3) X 

4 

(3 

X 2) 

X 

4 

6 

250 

185 

82 

PDS 364 

13:20:03.5 

-62:23:51.7 

2010-03-31 


300 

X 

4 


300 X 

4 


90 

X 

4 

3 

130 

116 

168 

PDS 69 

13:57:44.0 

-39:58:47.0 

2010-03-29 

(15 

X 2) 

X 

4 

(15 

X 2) X 

4 


5 

X 

4 

6 

45 

66 

149 

DG Cir 

15:03:23.4 

-63:22:57.2 

2010-03-31 


360 

X 

4 

(180 

X 2) X 

4 


10 

X 

4 

20 

14 

103 

138 

HD 132947 

15:04:56.2 

-63:07:50.0 

2010-03-12 

(20 

X 2) 

X 

4 

(20 

X 2) X 

4 


15 

X 

4 

3 

232 

265 

236 

HD 135344B 

15:15:48.2 

-37:09:16.7 

2010-03-31 

(5 

X 2) 

X 

4 

(5 

X 2) X 

4 


2 

X 

4 

6 

42 

116 

307 

HD 139614 

15:40:46.3 

-42:29:51.4 

2010-03-28 

(10 

X 2) 

X 

4 

(10 

X 2) X 

4 


5 

X 

4 

8 

35 

139 

127 

PDS 144S 

15:49:15.4 

-26:00:52.8 

2010-03-31 


300 

X 

4 


300 X 

4 

(5 

X 5) 

X 

4 

10 

49 

127 

124 

HD 141569 

15:49:57.8 

-03:55:18.6 

2010-03-28 

(2 

X 3) 

X 

4 

(2 

X 3) X 

4 

(2 

X 2) 

X 

4 

6 

177 

151 

144 

HD 141926 

15:54:21.5 

-55:19:41.3 

2010-03-12 

(15 

X 2) 

X 

4 

(15 

X 2) X 

4 


3 

X 

4 

10 

66 

198 

78 

HD 142666 

15:56:40.2 

-22:01:39.5 

2010-03-28 

(15 

X 3) 

X 

4 

(15 

X 3) X 

4 


2 

X 

4 

20 

53 

126 

114 

HD 142527 

15:56:41.8 

-42:19:21.0 

2010-04-01 

(5 

X 2) 

X 

4 

(5 

X 2) X 

4 


1 

X 

4 

20 

23 

88 

305 

HD 144432 

16:06:57.8 

-27:43:07.4 

2010-03-12 

(10 

X 3) 

X 

4 

(15 

X 3) X 

4 


2 

X 

4 

20 

55 

126 

109 

HD 144668 

16:08:34.0 

-39:06:19.4 

2010-03-30 

(3 

X 3) 

X 

4 

(3 

X 3) X 

4 

(0.75 

X 2) 

X 

4 

20 

96 

191 

276 

HD 145718 

16:13:11.4 

-22:29:08.3 

2010-03-29 

(15 

X 2) 

X 

4 

(15 

X 2) X 

4 


3 

X 

4 

10 

85 

165 

165 

PDS 415N 

16:18:37.4 

-24:05:22.0 

2010-03-31 


300 

X 

4 


300 X 

4 

(10 

X 5) 

X 

4 

5 

13 

66 

103 

HD 150193 

16:40:17.7 

-23:53:47.0 

2010-03-30 

(25 

X 3) 

X 

4 

(15 

X 3) X 

4 


2 

X 

4 

20 

100 

152 

322 

AK Sco 

16:54:45.0 

-36:53:17.1 

2009-10-05 

(10 

X 2) 

X 

4 

(10 

X 2) X 

4 


5 

X 

4 

8 

26 

91 

198 

PDS 431 

16:54:58.9 

-43:21:47.7 

2010-04-01 


300 

X 

4 


300 X 

4 


150 

X 

4 

2 

174 

171 

84 

KK Oph 

17:10:07.9 

-27:15:18.6 

2010-03-26 


200 

X 

4 

(100 

X 2) X 

4 

(2 

X 6) 

X 

4 

20 

58 

274 

231 

HD 163296 

17:56:21.4 

-21:57:21.7 

2009-10-05 

(3 

X 3) 

X 

4 

(3 

X 3) X 

4 

(1 

X 2) 

X 

4 

20 

87 

241 

127 

MWC 297 

18:27:39.7 

-03:49:53.1 

2009-10-06 


300 

X 

4 


10 X 

4 

(0.665 X 15) 

X 

4 

20 

162 

96 

127 


and 40 from the IVieira et al.l (l2000h catalogue, bring¬ 
ing the total number of targets to 91. The observations 
cover around 70% of the sou thern HAeBes identified by 
Ixhe. de Winter fc Perej lll99^. ari d about 50% of the tar- 
gets observed bv I Vieira et al. I^OOOl l. For many of these tar¬ 
gets little information is know n about them, particularly 
in regard to multiplicity (see IPuchenel l2015l . for a re- 
view). It is known t h at HAeBes have high binary fractions 
i jBahies et al.l l2006l : IWheelwright, Oudmaiier fc GoodwinI 
|2010D . and as a consequence, any close separation binaries 
will contribute towards observed spectra. In this work our fo¬ 
cus is on the UV and optical portions of the spectra, where 
we assume that the primary star, the HAeBe target, pro¬ 
vides the largest contribution to the brightness. Contribu¬ 
tions from secondary stars will be greater in the observed 
literature photometry as they use either larger slit widths 
or an aperture greater than the slit widths used here. Pho¬ 
tometry of all the targets are sourced from the literature, 
and are provided in Table [XI in Appendix 


Telluric standards were observed either just before or 
after each science exposure. They were observed in stare 
mode for short exposures, ~10 s, due to their brightness. 
Flux standard stars were observed on approximately half of 
the evenings in offset mode (offset mode allows accurate sky 
subtraction to be performed). 


2.2 Data Reduction 

All data were reduced foll owing standard proced ures of the 
X-Shooter pipeline vO.9.7 jModigliani et al1l2010li . Only one 
aspect is not included in the standard procedures and that 
is flux calibration. However, we do not require a flux calibra¬ 
tion for the work presented in this paper, we focus instead 
on determining the spectral shape. This spectral shape is 
needed in the UVB arm; specifically across the Balmer Jump 
region where the difference between the f7-band and J3-band 
is required in order to measure any excess UV-emission. The 
analysis of the spectral lines and the derivation of their lu¬ 
minosities is deferred to Paper II. 

Ordinarily, to obtain the correct spectral shape a flux 
calibration is performed using the observed flux standards 
of each night. However, as mentioned, flux standards were 
not observed on all evenings. A solution to this is to instead 
use the telluric standards, for which there is at least one per 
target, as a means of correction. This will ensure a uniform 
treatment to all of the targets. Caution should be noted of 
using non-flux standard stars for calibration; to mitigate any 
problems that could arise from this, consistency checks are 
made against the flux standards for the nights where they 
are available and will be discussed at the end of this sec¬ 
tion. For this method of spectral shape calibration accurate 
knowledge of the spectral type of each telluric is required. To 
ensure a homogeneous reduction we adopted our own spec- 
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Figure 1. Shown above are the two different cases of correcting target spectra via either: a flux standard - shown on the left; or a 
telluric standard - shown on the right. The top-left panel shows the observed flux standard (blue) and the true flux of the standard star 
(dashed-black). Similarly the top-right panel shows the observed telluric standard (green) and its expected spectra (dashed-black). The 
middle two panels show a division of the observed standard stars by their expected spectra from the top panels (blue and green). A two 
part curve (red) is also shown as a fit to this division. The bottom panel shows the result of applying the two fits from the middle panel 
to the target spectra (dashed-black). It can be seen that the two methods of correction are equivalent by lying on top of each other. All 
spectra have been arbitrarily scaled in order to be visible on each plot. 


tral typing of each telluric in this work. This helps to min¬ 
imise any reduction errors, and will also allow us to place a 
systematic error on this reduction method. Full details of the 
spectral typing, along with a discussion of how they compare 
with literature values, will be provided in Section [S] Once 
the spectral type is determined the observed telluric spec¬ 
trum is divided through by a model atmosphere of the same 
spectral type in order to obtain an instrumental response 
curve. The model atmospheres adopted here, and t hrough¬ 
out this work, are sets of Kur ucz-Castelli modelsjjKun^ 
Il993l: lUastelli &: Kurucd l2004l) computed by iMunari et al. 
( 2 nnRh . due to their small dispersion of 1 A over the UVB 
wavelength range (these will be referred to as KC-models 
hereafter). The resulting response curve from this division 


is then fit with two curves: one for the echelle orders where 
A < 3600 A and another for the orders where A > 3600 A. 
This is because the response at ~ 3600 A is not the same 
between the two over-lapping echelle orders. Figure [T] shows 
the procedure of the above method, for a target star for 
which both flux standard and telluric standards were ob¬ 
served, and highlights the two different response curves in¬ 
tersecting around the 17-band region in the middle panels. 
The figure also provides a consistency check by comparing 
the flux standard reduction, on the left, to the telluric stan¬ 
dard reduction, on the right. The bottom panel of the figure 
demonstrates the similarity of both results with a difference 
of < 3% across the spectra. Larger deviations are seen be¬ 
tween the two spectra close to 3000 A, due to low levels of 
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counts. This region is not used in this work and can be dis¬ 
regarded. This same check is performed on other stars for 
which both a telluric and flux standard are available, and the 
maximum deviation observed is only 5% across the spectra. 
Overall, it can be seen that this method of using the telluric 
for instrumental response correction provides a satisfactory 
calibration of the data, and is therefore performed on all 
targets. 


3 DETERMINING THE DISTANCE AND 

STELLAR PARAMETERS 

Determining accurate stellar parameters is crucial for ex¬ 
tracting an accretion rate, and for obtaining further infor¬ 
mation about the age, evolution, and ongoing processes in 
the environment around HAeBe stars. Many stars in this 
sample have had their stellar parameters determined previ¬ 
ously, but this has o ften been done in smaller sub-sets using a 
variety of me thods dMora et al.ll200ll:[ Iernandez_et_alj2004 
Mhjioi et al.l l2006l : iMontesinos et al. 20091 : lAlecian et al.l 
2011 see also the Appendix for additional references). For 
this reason a full treatment of determining stellar parame¬ 
ters is performed on the entire sample, in a homogeneous 
fashion, to provide better consistency between the stars. A 
comparison will also be made with the literature values to 
confirm the method employed; as most stars would be ex¬ 
pected to have similar temperatures to the previous litera¬ 
ture values. 

The determination of parameters is performed in a 
three step process: 1) Spectral typing is performed using the 
X-Shooter spectra to provide accurate limits on the effective 
surface temperature, Teff, and where possible the surface 
gravity too, log(g); 2) KC-models and the photometry are 
used to assess the reddening, Ay, and distance/radius, 
D/Ri,, ratio towards the targets; 3) Finally, PMS evolu¬ 
tionary tracks are used to infer a mass, M*, and age (and 
other parameters if not determined yet). The stages of this 
process are now given in detail. 


3.1 Temperature and Surface Gravity 
Determination 

The first stage takes advantage of the large wavelength cov¬ 
erage and good spectral resolution of X-Shooter to perform 
spectral typing; allowing us to narrow down the possible 
and log(g) of e ach target. This is done by following a sim¬ 
ilar method to IMontesinos et al.l (|2nn9h . of spectral typing 
using the wings of the hydrogen Balmer series, and the con¬ 
tinuum region 100-150A either side of the lines. These lines 
are favoured due to their sensitivity to changes in Tefr and 
log(g). Specifically, the H/3, Hy, and H5 lines are used as they 
have the largest intrinsic absorption of the series, except for 
the Ha line. Ha is not used for spectral typing as it is often 
seen entirely in emission, with the emission being both the 
strongest and broadest of the Balmer series in HAeBes. This 
could in-turn affect the derived parameters. Therefore, fit¬ 
ting of models to the line wings is performed using the other 
lines in the series. To perform the fitting, each line is first 
normalised based on the continuum either side of the line. 
They are then compared against a grid of KC-model spectra, 


which have also been normalised in the same way using the 
same regions either side of the line. The resolution of the grid 
is set to be in steps of 250 K for Tefr and 0.1 dex in log(g). 
The metallicity is kept at [M/H] = 0 throughout, although 
it has been shown that the choice of metallicity c an affect 
spectral typing in HAeBes (IMontesinos et aT1l2009ll . The fit 
of the synthetic spectra to the observed spectra is judged 
using the wings of each line, and continuum features, where 
the intensity is greater than 0.8; the line centre is excluded as 
it can often be found in emission. This approach avoids the 
problems of both emission, and rotational broadening in the 
line. Figure [2] gives four examples of this fitting, highlighting 
the power for obtaining an accurate Teg and log(g), where 
many errors are as small as the chosen step size. However, 
despite this reliable technique, issues arise for two cases. The 
first, is that there is a non-linear relationship between the 
Balmer line width an d the surface gravity fo r objects which 
have Tefl < 8000 K llGuimaraes et al.ll2()06lj . However, for 
temperatures up to 9000 K there is increased uncertainty 
due to the large presence of absorption features, which make 
normalising and comparing different surface gravity scenar¬ 
ios increasingly difficult. For these reasons we do not con¬ 
strain log(g) using the spectra for stars with a suspected 
Teff < 9000 K. 

The widths of the Balmer lines are tightly correlated to 
Teff and log(g), to the point where different combinations 
of the two can produce the same widths. However, this de¬ 
generacy can be broken when viewing the whole of the line 
profile and the absorption features within them (and also 
the photospheric absorption features outside of the wings). 

The second issue concerns objects which display very 
strong emission lines; where the line strength is exceptionally 
strong across the Balmer series to the point where the width 
of the lines eclipse even the broad photospheric absorption 
wings. Extremely strong P-Cygni, or inverse P-Cygni, pro¬ 
files can also affect the line shape in the wings. An example 
of extreme emission is shown in the bottom-right panel of 
Figure[2l where none of the intrinsic photospheric absorption 
lines can be seen due to the emission. P-Cygni absorption is 
also present in this example further complicating any pos¬ 
sible analysis of the wings. Objects, like the example just 
given, where both Teff and log(g) cannot be constrained by 
this method, will be treated separately on an individual ba¬ 
sis and are detailed in Appendix m The objects for which 
Teff has been constrained can have all of their parameters 
determined in the next two steps. 

For the telluric standards the same above steps are ap¬ 
plied. This is because they are well-behaved stars for which 
a Teff and log(g) determination is straight forward. These 
parameters are required for the data reduction discussed 
previously in Section [2.21 

Figure |4] compares the temperatures derived in this 
work against previous estimates from the literature (see Ta¬ 
ble |A] in Appendix The temperature is chosen for com¬ 
parison as it is a key stellar parameter which can be deter¬ 
mined more readily than log(g), and its appearance in the 
literature is more frequent than other parameters (allowing 
a greater number of comparisons to be made). The major¬ 
ity of literature works provide a spectral type rather than a 
precise temperature so we assign an error of 10% for these. 
The figure shows over 95% of the stars are in agreement, 
within the errors. Also, the temperature determinations in 
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Table 2. Details of the derived stellar parameters. Stars for which a distance estimate from the literature is used, or if the 
star was moved to the ZAMS, are noted in the final column (Column 9, a legend is provided as a footnote). 


Name 

TeS 

log(g) 

log(i*) 

M* 

R* 


Age 

Distance 

Notes 


(K) 

[cm/s^] 

[L©] 

(Mq) 

(R©) 

(mag) 

(Myr) 

(pc) 


UX Ori 

8500 ± 250 

Q QQ+0-25 
o.yu_Q 25 

1 54+0-38 
3.04_q 33 

2 - 1 +i 

9 7 + 0.4 
' - 0.2 

n 48+0-07 

d-4o_o.o3 

4.24lljt 

600l°° 


PDS 174 

17000 ± 2000 


2 

"^•^-^-0.56 

5.o1|3 

q q+0.7 

o.o_o 7 

q Cl+0.07 
'3-O1-0.07 

0.60l°;« 

11261238 


V1012 Ori 

8500 ± 250 

4.38l°;l| 


l-'0-0.3 

1 4+0-4 
1 -^- 0 .4 

1 q 9 + 0.02 
i.OZ_Q Q 4 

15.161^:1 

445+137 
^40 _ 1^30 

t * 

HD 34282 

9500 ± 250 

4 40"*"*^'^^ 

1 1 7+0.28 
^•-^'-0.36 

1 9+0-4 
1 -^- 0 .4 

1 4+0.4 
1 -^- 0 .4 

0 01+0 02 
o'-o'i-o.oo 

10 . 00 l |;00 

360+111 

'300-109 


HD 287823 

8375 ± 125 

4 

l.09l°;i^ 

1 -7+0.1 

i-'-o.i 

1 7+0.3 
l-'-0.3 

6 '- 6 IU_o 00 

9.0112.28 

340l«® 

t a 

HD 287841 

7750 ± 250 

4 

^•^'- 0.12 

0.87l°;i^ 

-*^•^- 0.1 

1 5 + 0.2 
1 -^- 0.2 

6 '- 6 IU_o 00 

14 07+4-10 

1^-9(_4 10 

340l°® 

t a 

HD 290409 

9750 ± 500 

^•^^-0.25 

1 42+0-29 

9 1 +0.2 
^• 1 - 0.2 

1 84"0.1 
^•^- 0.1 

6 '- 6 IU_o 00 

5 50+2.03 
O.OU _2 03 

5141°® 


HD 35929 

7000 ± 250 

3.47l°;“i 

1 70+0-17 

2 9+0-4 
^•^-0.4 

5 2 + 0.7 
^•^- 0.8 

0 nn+^’^^ 
U.UU_o.oo 

1 03+0-81 

I-OD-0.55 

36 OI 72 

tb 

HD 290500 

9500 ± 500 

3.80l°i° 

1 Q4+0-01 

l-y4_0.53 

9 8+1'^ 
^•^-0.7 

0 5 + 1.0 
'0-^-0.5 

0 00 +*^'*^^ 
U.UU_o.oo 

2.26^172 

1522l«6 


HD 244314 

8500 ± 250 


1 91 +0.18 
-*^•^-*^-0.23 

1 R+^’l 
-*-•^- 0.1 

1 9 + 0.3 
l-l^-0.3 

0 1 0 + 0.02 
U.lU-0.05 

7 59+1-97 
< •9^_1 96 

440li 

t a 

HK Ori 

8500 ± 500 


1 1 ‘3+0.24 
-L-J- 0 _q 27 

1 7 + 0.2 

1 - ‘ - 0.2 

1 7 + 0.3 
l-'-0.3 

1 91 + 0.12 

l-^l-0.14 

8 70 + 2.73 

0 . ' 0_2 73 

440l«« 

t a 

HD 244604 

9000 ± 250 

o.y»_Q 13 

1 54+0-18 
l-94i_Q 23 

9 1 +0.3 
^• 1 - 0.2 

2 4+0-^ 

^•^-0.5 

0 14+0-00 

4.56tl:g 

440l«« 

t a 

UY Ori 

9750 ± 250 

^.ou_Q 20 

1 30+0-23 

i.ou_o 23 

2 0 ^^’^ 
^•'^-O.l 

1 7 + 0-1 
i-'-o.i 

111 + 0.02 
1 - 11 - 0.00 

6.35+il 

10271°® 


HD 245185 

10000 ± 500 

^•^^-0.25 

1 49+0-29 

-L-^»_0.29 

9 9 + 0.2 
^•^- 0.2 

1 84"0.1 
^•^- 0.1 

0 00 +^’^^ 

4 91+^'^0 
^•^1-1.70 

5191®" 


T Ori 

9000 ± 500 

Q «n+0-30 
o.DU_o 30 

9 1 9+0.47 

q q+1.5 
'^•'^- 1.0 

A 84"1.0 

4-^-0 .8 

l-OU_o.o5 

1 qc:;+2.27 
J^.OO_o.93 

750459 


V380 Ori 

9750 ± 750 

^•UU_o.35 

1 71+0.52 

1 -' 1-0.26 

2 3+11 
z.O- 0.2 

9 c;+0.5 
^• 0 - 0.1 

9 91 +0.05 
"^•^1-0.07 

q 7q+2.46 

0 . <o_2.49 

330 I 1 ? 


HD 37258 

9750 ± 500 

4.ZO_o_25 

1 42+^-^^ 

l-4-i_0.29 

9 1 + 0.2 
^• 1 - 0.2 

1 84"0.1 
1 -^- 0.1 

0 06+*^'*^^ 
U.UD_q Q 4 

5 50+2.03 
O.OU -2 03 

424l°5 


HD 290770 

10500 ± 250 


1-9^_0.38 

9 q+0.5 
^•'1-0.5 

^•O'- 0.2 

6 '-»JU_o 00 

^-10-1.89 

4401^5 


BF Ori 

9000 ± 250 

Q Q7+0-15 
'J-y < -0.13 

1 r7+0.19 

l-^'- 0.22 

9 1 +0.3 
^• 1 - 0.2 

9 ^+0.4 

^•^-0.5 

r\ qq+0.03 

u.oo _0 02 

4.34+:g 

51 0+102 
01 U_io 2 

t a 

HD 37357 

9500 ± 250 


1 r 9 + 0.13 

2 1+^1 
^• 1 - 0.1 

2 1 + 0-1 
^• 1 - 0.0 

6 '-»JU_o 00 

4 93+0-87 
^.»o_o S7 

344+15 


HD 290764 

7875 ± 375 


i.90_Q 26 

1 9+0-4 
1 -^- 0.2 

^•'0-0.5 

1^-19-0.14 

c 9 c: + 2.58 

O-^O-i 90 

470l°l 

t a 

HD 37411 

9750 ± 250 

^•oO_o.i5 

1 9Q+0.28 

i.zo_Q 36 

1 9+0-4 
1 -^- 0 .4 

1 5 + 0.5 
1 -^- 0 .5 

0 21 + 0-01 

9.^ 1 - 0.00 

g.ootljo 

qr8+109 

'300-107 


V599 Ori 

8000 ± 250 

o ■tq+0.13 
'^- 0.12 

l- 68 lgil 

9 5 + 0.4 
^•0-0.4 

3 «+ 0.6 

O.D_o 7 

4 05+0.06 
^•00_Q Q7 

9 09 + 1.40 
^• 0 ^- 1.00 

51 0+102 
OlU -102 

t a 

V350 Ori 

9000 ± 250 


1 oi +0.19 
i-'^i- 0.22 

1 9+0-1 
1 -^- 0.1 

1 9 + 0.3 
l-l^-0.3 

0 09+0.02 
U.Dy_Q Q3 

0 41 +1-97 
'^•41-1.67 

51 0+102 
OlU -102 

t a 

HD 250550 

11000 ± 500 

3-80l°i° 

9 90 + O. 6 I 
z.^o_o 53 

3 4 + 2.1 

q 0 + 1.0 

'=^•^-0.5 

0 00 +*^'*^^ 
9-UU_o 00 

1 42"*”^'^^ 
i-^^-1.09 

973+267 

»<'J-136 


V791 Mon 

15000 ± 1500 

16 

9 05 + 0.40 

^•90_q 51 

9.0_q 7 

9 9 + 0.7 
^•^-0.7 

1 1 7+0.06 
l-l'-0.04 

1 80+0-90 

l.o^_0.90 

648l?°l 


PDS 124 

10250 ± 250 

^.ou_o 20 

1 47+0-23 
1-^' -0.23 

9 9 + 0.2 
^•^- 0.2 

1 7 + 0.1 

i-'-o.i 

1 9q+0.03 

1.^0_Q QQ 

c- 40+1.64 

5-40-i 64 

894l®4 


LkHa 339 

10500 ± 250 

^•^^- 0.20 

l-^^-0.24 

9 q+ 0.2 
^•' 1 - 0.2 

2 .ot°:l 

3 54+0-01 
9-04i-o 01 

4.i6+:l^ 

5971°^ 


VY Mon 

12000 ± 4000 

o yir+0.50 

O- i 'J_0.50 

z.oo_o.65 


4.4+:? 

5 08+0-17 

O.DO_o.45 

0.89lg;5J 

43911®! 


R Mon 

12000 ± 2000 

4.UU_q 24 

2 19+0-34 

-i-iy_0.40 

q 1 + 0.8 

'^•i-O.O 

9 q+0.3 

2 42+0.08 
^•^^-0.13 

1 92+^'^^ 
i-y^-0.89 

80011°° 

t c 

V590 Mon 

12500 ± 1000 


^•OID-0.37 

q 1 + 0.6 
'^• 1 - 0.6 

9 q+ 0.2 

Z.O_o _2 

1 03+0.04 

l-UO_o.o5 

2 19+1-00 
^• 1 ^- 1.00 

1722ll^J 


PDS 24 

10500 ± 500 

^•^^-0.30 

i-9^_0.38 

9 q+0.5 
^•'^-0.5 

2 0 + 0.2 
^•O'- 0.2 

1 11+0.04 
1-11-0.03 


164611®? 


PDS 130 

10500 ± 250 

o.»u_o 20 

9 r »9 + 0.30 
^•'3^-0.27 

9 8+0-7 
^•^-0.4 

0 1 + 0.4 

'O- 1 - 0.2 

2 07+0-01 
- 0.00 

9 9c;+1-32 

1748111® 


PDS 229N 

12500 ± 250 

^•^'^- 0.20 

Z.UO_o 23 

q 1 + 0.2 
'^• 1 - 0.2 

9 q+ 0.1 

^•'O-o.l 

2 0 O+O.O 1 

Z.UO_Q QJ 

2 19+0-63 

^•1^-0.63 

1379151 


GU CMa 

22500 ± 1500 

o.yu_o 40 

q 07 + 0.62 
'3-o< _o.62 

9 5 + 6-8 
^• 0 - 6.8 

5 7 + 1.7 

0+-1 7 

0 57+0.03 

9-0^-o 02 

0 . 1 +rio 

cqi +163 
0'^i-163 


HT CMa 

10500 ± 500 

4.UU_q 20 

1 qq + 0.29 

i.05_Q 24 

2 . 6 t°:® 

2 0+0.3 

^•'O-O.l 

r\ 9q+0.05 

U.ZO_Q Q2 

2.96tii^ 

163411®'! 


Z CMa 

8500 ± 500 

9 

z.oo_Q_i7 

o.uz_Q 28 

ii-o+i:? 

29.815.2 

q q 7 + 0.12 

9.0(-o.i6 

0 0O+O.O2 
U.UO_Q Q2 

1050l?l° 

td 

HU CMa 

13000 ± 250 

^•^^- 0.20 

2 10 "^*^'^^ 
-i-J-U_Q 23 

q 9 + 0.2 
'^•"^- 0.2 

2 4 + 0.1 
^•^- 0.1 

0.80l°;°? 

1 88+0-52 
l-oo_Q_52 

12404 ! 


HD 53367 

29500 ± 1000 

^•^^-0.25 

4 11 + 0.37 
^•11-0.37 

12-3+i 

4.3lgl 

1 88 + 0-02 
l-oO-O.Ol 

o.08l°;°| 

340l|l 


PDS 241 

26000 ± 1500 


4 1 1+0.46 
^•11-0.46 

ii-e+j 

5.6li;2 

^-9U-o 01 

0.08t°:°i 

29071®!? 


NX Pup 

7000 ± 250 

o 1^0+0.13 

1 90 + 0.20 
-*^•^^- 0.21 

I'Slaa 

3 0+0-3 
'0-0'_o.5 

9-UU_o 00 

4.9212.37 

410li 

t a 

PDS 27 

17500 ± 3500 

3.16t°oil 

4 39+0.40 
^.o»_o 40 

15.3111 

17 0+4.0 

l'-0'-4.0 

5 0O+O.13 

9-U9_o 13 

0 10 + 0-10 
' 0 -i' 0 _o.io 

31701®?° 

t e 

PDS 133 

14000 ± 2000 

4.08l°;l? 

^•40_o.38 

q 7 + 0.9 
'^•'-0.7 

2 9 + 0.4 

1 43+0-09 

l-4o_o.io 

1 97+0.94 

l-^'-0.56 

25001®°° 

tf 

HD 59319 

12500 ± 500 

o.ou_Q 20 

o.uo_Q 3 Q 

5 7 + 1-6 
OW-i 2 

7 0+0.9 
' • 0 '- 0.8 

0 00+0-05 

u.uu_o.oo 

o.32l°;?^ 

121811 ®? 


PDS 134 

14000 ± 500 

o.4iU_o 30 

3 45+0-46 
9.40_q 45 

7 «+3.5 
'•'0-2.3 

9 1 + 1.0 
^•1-1.5 

1 99 + 0.03 

J^-^-i_ 0.02 

o.i5lg:?§ 

56871“!® 


HD 68695 

9250 ± 250 

4 40"*"*^'^^ 

111 +0.28 
1-11-0.37 

1 8+6-4 
^■®-0.4 

1 4+°-^ 
1-4^-0.4 

9-UU_o 00 

10.0015;°° 

34411?? 
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Table 2 — continued 


Name 

Teff 

log(g) 

log(D*) 

M* 

R* 

Av 

Age 

Distance 

Notes 


(K) 

[cm/s2] 

[Lq] 

(M©) 

(R©) 

(mag) 

(Myr) 

(pc) 


HD 72106 

8750 ± 250 

3.89l°;« 

1 fi3+0-18 

i.DO_o,21 

9 q + 0.3 
^•'^-0.3 

9 O+0.5 
^-^-0.5 

0 00+0-05 

U.UU-o 00 

3 76+^''^^ 
'^•<0-1.18 

370111 

t a 

TYC 8581-2002-1 

9750 ± 250 

4.00«;1° 

1 71+0-14 
' -'--0.12 

9 q+0.2 
^•'^-0.1 

2 5+°-l 

0 94+0-04 
^•^**-0.00 

q 7q+0.78 

•J- <0-0 89 

9021!! 


PDS 33 

9750 ± 250 

4.40l°;l| 

1 Qq+0.27 
i.zo_0 36 

1 Q+0-4 
1-^-0.4 

1 4+°-“^ 
1-^-0.4 

0 C9+0.04 
U.QZ-o 00 

9 nn+^'^0 

y.UU_4 50 

Qq 9 + 282 
yoZ-279 

* 

HD 76534 

19000 ± 500 

^•-‘-'-'-0.20 

q -| Q+0.26 
'^•-'^^-0.20 

6-olS:^ 

3 6+0-8 
'^-0-0.2 

0-621H? 

0 37+0.19 

O-O<-0 14 

568l« 


PDS 281 

16000 ± 1500 

O.JU_Q 3Q 

'^•0^-0.45 

8.3ll° 


1.89ino 

0 1 9+0.23 
O-J-^-0.09 

936ll°! 


PDS 286 

30000 ± 3000 

4 

^•^O_0.16 

4 10+0.41 

4.IO-0.52 

13.5^1^ 


c 97+0.05 

O.^/_0.04 

0 10+0.05 

U.lU-o 05 

521IIS 


PDS 297 

10750 ± 250 

4.UU_q 20 

1 03+0.29 
i.yo_o.24 

2 6+0-0 
^•0-0.2 

9 ^7+0.3 
^-'-0.1 

0.8ll°;°l 

9 77+1.12 
'-1.22 

14651!®® 


HD 85567 

13000 ± 500 

'^•^^-0.30 

q 1 q+0.46 
'^•^‘^-0.45 

0-^-1.8 

7 9+1-5 
‘ -^-1.2 

0-89l2:ol 

n 97+0.52 

O-'^<-0.18 

907111! 


HD 87403 

10000 ± 250 

o.ou_o 10 

9 Qq+0.15 
^•^'^-0.15 

5-+0.6 

8-+0.5 

0 00'*'0 08 
'O-l'O-O.OO 

0 32+0-16 

o.o^-O 10 

180111!! 


PDS 37 

17500 ± 3500 

^•^^-0.35 

4 73+0.39 
^O_0.39 


25.8l®;° 

c 01 +0.13 
O-ol-o 13 

0 lO+O-lO 

o-iO-O.io 

4310+070 

2±oiU-67o 

t e 

HD 305298 

34000 ± 1000 

4 O1+0.16 

^•^o_o.4i 

15.7111 

4.61M 

±.ou-o 02 

0 02 + 0-01 
o-o^-O.oi 

33661“!® 

* 

HD 94509 

11500 ± 1000 

2 90'''^'^^ 
^•^'^-0.40 

q 7C+0.65 
"J- '^-0.62 

10 . 814-3 

19. 2 I®;® 

U.UU-O 00 

0 0c:+O.16 
O.UO_o 05 

4384+^085 

^00^-1009 


HD 95881 

10000 ± 250 

Q on+O-lO 

2.98l°;lt 

6 . 2 l 0 ;« 


0 00+0-05 

U.UU-O.00 

0 21 +0.10 
O-^l-0.07 

1290 II! 


HD 96042 

25500 ± 1500 

Q or»+0.20 

4.oU_q 20 

4.36l°J| 


7 q + 1-3 
'-^-0.8 

n 70+0.03 

U. 'O-o.Ol 

o.02t°;“® 

1792 II!! 


HD 97048 

10500 ± 500 

4 ■^n+0.20 

4.c>U_o.20 

1 K9 + O .23 
i.OZ-0.23 

9 9 + 0.2 

Z.Z-o 2 

1 7+0.1 
l-'-O.l 

0 Q0+O-O5 

u.yU-0.02 


171 II 


HD 98922 

10500 ± 250 


^•"10-0.15 

4 0+0.5 
^•^-0.5 

K 9+0.3 

O-^-0.3 

U.UU-O 00 

0-84tg;l 

3461!! 


HD 100453 

7250 ± 250 

4.08t°;« 

0 93+0-17 

1 c+0.2 

I'O-O.l 

1 q+0.3 
1 -^- 0 .3 

u.uu_o 00 

9 97+0-20 
y.y / _2.79 

I 22 I!! 

fb 

HD 100546 

9750 ± 500 

4 34+0.06 

-'-•'^^-0.14 

1 Q+O-l 

I'^-O.l 

1 - 0 - 0.1 

u.uu_o 00 

7 02+^'‘^0 

<•0^-1.49 

9711! 

fb 

HD 101412 

9750 ± 250 

^.ou_o 20 

1 36+0-23 

±.ou_o 23 

^•^-0.1 

1 7+0.1 
l-'-O.l 

0-2llo.oo 

6 33+1-84 
o.oo_i g 4 

301 II! 


PDS 344 

15250 ± 500 

4.30l°i° 

9 qQ+0.25 

Z.c>y_o.25 

q 7+0.5 
'^•^-0.5 

9 q+0.1 

^-O-O.l 

0.86l°;°l 

1 4.8+0-62 
1-4O-0 52 

275611!! 


HD 104237 

8000 ± 250 

3.89l°/J 

1 41+0.17 
1-41-0.21 

9 q+0.3 
'^•^- 0.2 

^-0-0.4 

0 O0+O-O5 

u.uu_o.oo 

4 92+1-0’^ 

4.y^-i 46 

115l!| 

fb 

V1028 Cen 

14000 ± 500 

q Qn + 0.30 
-:>-oU_o.3o 

9 oc+0.45 

^•00-0.41 

4 7+2.0 
^•'-1.0 

4 r+0.9 

^-0-0.5 

0 1^7+0.01 
^•u<-0.03 

0 ko+0.69 

O-OU-o 39 

18431!®! 


PDS 361S 

18500 ± 1000 

3-80l°i° 

3 e:3+0.46 

'J-^'J-o 41 

7 4+8.2 
' '^-1.7 

K 7+1.1 

O-'-0.7 

i.yU-o 01 

0 19+0-23 

O-iy-0.12 

4385111! 


HD 114981 

16000 ± 500 

O.OU_Q 20 

3 47+0.30 
-0.30 

7 q+2.0 

'•'^-1.5 

7 1 +0.9 
'-1-0.8 

y-yu-o.oo 

n 1 q+0.17 

b.lO-o 09 

9081!“ 


PDS 364 

12500 ± 1000 

^•^'^- 0.20 

2-UO_o 23 

q 1 + 0.2 
'^•-'^- 0.2 

9 q+0.1 
^-■^-0.1 

1 07+0.05 
l-o<-0.03 

2-19ll!®58 

17151!! 


PDS 69 

15000 ± 2000 

^•^'^-0.35 

9 79+0.52 

2- ‘ ^-0.76 

4 3+2.0 

3 4+0.7 
'^-^-0.6 

1.0U_Q Q7 

0.841202 

63011" 


DG Cir 

11000 ± 3000 

4 41+0-18 

1 49+0-68 
l-^y_0.93 

9 9+0.4 

^•^-0.4 

r:+0.5 
1-0-0.5 

8 94+0-18 
>J-y^-0.54 

3 93+2-97 
o.yO-2.97 

71 q+250 
< J^'J_184 

t * 

HD 132947 

10250 ± 250 

3 QQ+O-IO 

o.yu_Q 10 

1.97l°;l® 

9 »7+0.3 

^•'-0.3 

q 1 + 0.2 
0 - 1 - 0.1 

0 O0+O-O5 

U.UU-O.00 

2 44+0-77 

5651®! 


HD 135344B 

6375 ± 125 

3 Q4+0-12 

0.85l°i» 

1 c:+ 0.2 

1-O-0.2 

9 9+0.4 

^-^-0.4 

n 9q+0.05 

u.zO-0.06 

7 99+8.24 
< •yy-2.34 

140l!! 

t a 

HD 139614 

7750 ± 250 

4 31+0-12 

0.82l°i7 

I'O-O.l 

1 4'*"0-^ 

1-^-0.2 

u.uu_o 00 

15.64lt;^® 

1401!! 

t a 

PDS 144S 

7750 ± 250 

4 13+0.14 

i-^^-0.23 

i-O-O.l 

1 q+0.3 
^-^-0.3 

0 C7+0.07 

y-o/-o.o8 

9.45ll,« 

loool!!! 

ff 

HD 141569 

9750 ± 250 

^•'50_o 15 

1 90+0.28 

1-4O-0 37 

1 0+0.4 
^•^-0.4 

1-+0.5 

0 Ol+O-Ol 

y-yi-0.00 

9 OO+^'OO 

y-OO-4.50 

1121!! 


HD 141926 

28000 ± 1500 

q 7C+O.25 

'O_o.25 

4 70+0-^6 
‘^- 0.37 

iy-4-5.0 

9-+1:3 

^-^01-0.04 

0.00l°;“ 

1254111! 


HD 142666 

7500 ± 250 

4 13+0.11 
^•1*^-0.16 

0-96to.24 

1 - 0 - 0.1 

1 q+0.3 
1-^-0.3 

0 C 0 +O.O 8 
U.OU_Q 09 

10.431®;=! 

1451!! 

t a 

HD 142527 

6500 ± 250 

q Qq+0.08 

o.yo_0.08 

0-90tg;li 

1 - 0 - 0.1 

9 9+0.1 
-^-^- 0.2 

0 O0+O-O5 

y-uu-o.oo 

8.0811;®! 

140l!® 

t g 

HD 144432 

7500 ± 250 

4.UO_0.i4 

l-O^-0.21 

1 6 + 0-2 
1 - 0 - 0.1 

2 0+0-8 
^•0-0.3 

y-ou-o.oo 

Q 79+4.81 

'^"^-2.50 

1601!! 

tb 

HD 144668 

8500 ± 250 

q 7 i::+0.13 
'J- <O _0 12 

1 76+0-19 

1- '’^-0.22 

9 c:+0.4 
^•0-0.4 

q c:+0.6 

O-O-0.6 

r» qq+0.05 
U.OO-O 04 

2-7011;!! 

1601!! 

tb 

HD 145718 

8000 ± 250 

4 3'7+0.15 

^•'^'-0.15 

n 09+0.29 
9-o2-o 37 

1 r+0.3 
1-0-0.3 

1 3+0.4 

1-O-0.4 

0 74+0-06 
'-'•'^-0.05 

19.54l!;ll 

1341S 

t * 

PDS 415N 

6250 ± 250 

4 47+0.15 

^•^^-0.15 

0 13+0-30 

1 1+0-2 
1-1-0.2 

1 0+0.3 
l-l'-0.3 

1 11+0-11 
1-11-0.15 

33b 02+108.01 

OOD.UZ_j^gg 04 

1971!! 

t * 

HD 150193 

9000 ± 250 

^•^'-0.17 

1 91 +0.19 
^•^^-0.23 

1 Q+0-1 

i-i'-o.i 

1 7+0.3 
1-^-0.3 

1 cc+0.02 

1-00_Q 04 

7 99+1.89 

' -^"^-1.89 

1201!1 

tb 

AK Sco 

6250 ± 250 

^•^o_o 10 

0.38l°l® 

1-^-0.1 

1 q+0.2 

1-O-0.2 

0 00+0-05 
U.UU_0 00 

17 71 +‘l-7l 
!'• ^1-3.42 

1031!! 

tb 

PDS 431 

10500 ± 500 

q 7/1+0.20 
''-'-0.20 

9 q9 + 0.31 

Z.OZ_Q 3Q 

3 c+l.O 
9.0-07 

4-4tS:5 

1 76+0-03 

1- <0-0 03 

i-i9lS;!I 

28751®!! 


KK Oph 

8500 ± 500 

4.38l°;l| 

0 Q4+0-33 

U.y4_0 43 

1 q+0.3 
1-0-0.3 

1.4tH 

2 70+0.10 
^•<'1-0.15 

I5.i6ll;«« 

279+!! 

t * 

HD 163296 

9250 ± 250 

4.30t°;2° 

1 9q + 0.23 
x.zo_Q 23 

1 Q+0-1 
1-^-0.1 


0 OO+O-O^ 

U.UU-O 00 

7.561111 

loii! 


MWC 297 

24500 ± 1500 


3 03+0.46 
o.yu_Q 4 g 

lO-^-4.6 

5.311:1 

0 47+0.04 
o-‘^<-0.03 

o-iolH! 

i7oi!! 



t - A literature distance is initially adopted to th ese stars, as 1 ok ( s) can not b e determined f rom t he sp ectra alone. * - Stars w hich 
have been placed on th e Z AMS. Reference s : fa^ Ide Zeeuw et al, 

ISheychenko et al\ l|l999h . (e) lAbabakr et alJ l|2QlV 


_ Il999^. (b) Ivan Leeuwei ] ll2007t ). (cj [ Pahm Sz Simonl ll2005lh I'd') 

accepted), (f) IVieira et al.l^Q03h . fg'i iFukagawa et alJ l|2006l 'l. ('h'l iLoinard et alJ 
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this work have been based on some of the best spectra avail¬ 
able for these objects, which helps keep errors to a minimum. 
This serves as a justification for the homogeneous approach 
to determining temperatures and their use here, for both the 
target stars and the telluric standards alike. 


3.2 Photometry Fitting 


The second step of this process takes two directions: One 
case is where both Tefr and log(g) could be determined from 
the spectra, and the other case is for when only Tefj could 
be determined. In both cases fitting spectra of model atmo¬ 
spheres, based on the parameters determined in the previous 
step, to the observed optical photometry will be performed. 
The fitting will provide a level of reddening, Av to each star, 
and a scaling factor, D/Ri,, due to the fitting of surface flux 
models to observed photometry. An accurate temperature is 
paramount here in order to break any degeneracy of fitting 
models to the photometry. 

To perform the fitting only the BVRI points are used; 
the [/-band can be influenced heavily by the Balmer Ex¬ 
cess, and no photometry long-wards of the 7-band is used 
due to the possible influence of the IR-excess (which itself 
would require dedicated modelling). The 73-band can also 
be affected in the cases of extremely large flux excess. For¬ 
tunately, these cases are rare and the change in the i3-band 
magnitude would not significantly affect the fitting (the fit¬ 
ting is far more sensitive to the input temperature). 

Another point to consider when looking at opti¬ 
cal photometry is the effects of variabil ity, as this has 
been observed in numerous HAeBes (|de_WinteL_et_aD 


2001 I: Oudmaiier et al.l 2001: _ Mendieutfa et al. 2011a 


Pogodin et al.ll2012l : lMendigutfa et al.ll20^13li . However 


ability information is not present for all of the targets, but 
we estimate that the calculated parameters will not be af¬ 
fected significantly if the photometric variation is less than 
0.2 mag. In all cases we use photometry when at maximum 
brightness, as this best reflects the scenario where we are 
mostly viewing the stellar photosphere. So an assumption is 
adopted here that the photometry we use is predominately 
photospheric and not highly variable. 

In order to fit the photometry a unique grid of KC- 
models is set up based on the limits derived in step 1 for 
each star; the grid follows the same step sizes used in the 
previous step too. Log(g) does not have a significant effect 
on the fitting to the photometry, as the spectral shape is 
overwhelmingly dominated by the temperature. This allows 
log(g)=4.0 to be adopted and used in this step for the stars 
where log(g) could not be determined from the spectra; this 
value will be revised in the next step. 

The models are reddened until a best fit to the pho¬ 
tometry is achieved; the best fit being when the red¬ 
dened SED shape of the model is in-line with the pho¬ 
tometry. The dereddening is performed usin g the redden¬ 
ing law of ICardelli. Clayton fc Mathid (Il989l l. with a stan¬ 
dard Rv=3.1, in all cases. It is possible that the total-to- 
selective extinction may be higher, possibly Rv = 5.0, for 
some of th e stars in this sample based on previous a nalysis 
of HAeBes jHernandez et al.l2004l : [Manoi et al.l2006l l. How¬ 
ever, the choice of Rv will only affect the targest with the 
most extinction and the changes this will have on the stel¬ 
lar parameters and accretion rates are minimal due to the 


observed colour excess remaining the same. The majority of 
the targets have a mean 75(73 — F) « 0.4. 

Returning to the fitting, the model is normalised to the 
F-band point by a scaling factor, which is (77/73*)^. This 
scaling factor arises from the fitting of models in units of sur¬ 
face flux to observed photometry. An advantage of knowing 
this scaling factor is that it allows either distance or radius to 
be determined provided the other is known. Figure [3] shows 
an example of the above fitting for the case of V1012 Ori, 
along with a dereddened version of the photometry and the 
model spectra. This object is shown as it demonstrates a 
clear IR-excess, a noticeable Av, and a [/-band magnitude 
slightly higher than the KC-model spectra (possible Balmer 
Excess). 

At this stage the techniques diverge between the stars 
for which a log(g) was determined, and for the ones in which 
it could not be. For the former no further action is taken in 
this step. For the latter a distance is adopted to the star 
based upon the location of the star on the sky and its possi¬ 
ble associations with nearby star forming regions. The stars 
for which this is performed are noted in Table [21 the litera¬ 
ture distances adopted and references are both provided in 
the same table (and also in Table in Appendix O . An 
error of 20% is adopted for the distance, as this helps reflect 
the additional uncertainty on whether the star is truly part 
of the association, and the possible extent of the associa¬ 
tion. If the error is higher than 20% then the higher error 
is adopted instead. By adopting a distance to these stars 
a radius can be determined from the scaling factor. Then, 
combining this radius with the temperature, the luminosity 
is calculated by a black-body relationship of L* = 

(this calculation is equivalent to the sum of the flux under 
the KC-model multiplied by 47r77^). 


3.3 Mass, Age, Radius, and log(g) Determination 


In this third and final step, the remaining stellar param¬ 
eters are now determine d through the use of PMS tracks. 
The PARSEC tracks of iBressati et al.l 1I2OI2I') are used for 
the majority of this step as they cover a mass range of 0.1- 
I 2 M 0 , which encompasses all of the theoretical HAeBe mass 
range, and a metall icity is chosen of Z =0.01 (this is close 
to solar met allicity llCaff au et al .]l201llH. Additionally, two 
tracks from iBernasconi fc Maeder (Il996l') are used for ob¬ 


jects greater than 12 Mq. Each track is of a fixed mass, 
with no accretion contribution, which evolves over time in 
Tefi and L* as the star contracts. As TeS and L* change so 
do 73* and log(g) as a consequence. This allows each star 
to be plotted on either: a L* vs. Teff set of tracks, for the 
stars where L* is known from the adopted distance; or on a 
log(g) vs. Tefj set of tracks, for the stars where both log(g) 
and Tefi were determined from the spectra. For the first sce¬ 
nario a mass and age are extracted from the PMS tracks. 
Then, log(g) is calculated using this mass and the radius 
from the previous step. For the second scenario, luminos¬ 
ity, mass, and age are all extracted from the tracks. These 
can then be used to obtain a radius from the temperature 
and luminosity; or the mass and log(g), both choices are 
equivalent. Finally, using the 77/73* factor a distance can be 
determined. 

However, not all cases allow parameters to be extracted 
from the tracks. These few cases are where the stars are lo- 
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Wovelength (A) Wavelength (A) 

Figure 2. Examples of spectral typing for four targets are presented here. Each panel shows both the target spectra (grey), and a 
KC-model which denotes a good fit (red). The parameters for the KC-model are given for each fit. Also plotted is a dashed line, at 0.8 
of the normalised intensity, which is used as a cut-off in the fitting. The two panels on the left show the cases for straight forward fit 
where there is no obvious emission. The figure in the top-right is a case where there is clear emission present; a good fit is still achieved. 
The bottom-right panel gives an example of one of the few objects which cannot be spectrally typed in this way due to extremely strong 
emission; this exceptional object, and others like it, are presented in detail in Appendix iBl 


Gated below the zero-age main sequence, ZAMS. It shonld 
be noted that for a few of these cases, where the stars are 
only just below the ZAMS of the chosen tracks, then tracks 
relating to stars with a lower metallicity may be more ap¬ 
propriate. However, in general, it appears more likely that 
their placement is genuinely below the ZAMS and is dne 
to the adopted literature distances used being incorrect; as 
their use provides small radii from the I?/W* ratio. The ra¬ 
dius is deemed too small as it is less than the expected ra¬ 
dius of a ZAMS star of the same temperature. Addition¬ 
ally, most of these stars have Diffuse Interstellar Bands, 
DIBs, in their spectra whi ch suggest Ay ~ 0.5 — 2.0 mag 
djenniskens fc Desert|[l9S^ . Extinction due to DIBs follows 
a trend of ~ 1.8 mag/kpc dWhittetl l2003h . This suggests 
that the distances should be greater than the adopted val¬ 
ues and should be revised. Previously, for these stars, the 
assumption had been made that the stars are associated 
with a star-forming region. It is now more probable from 
the spectral typing and position of the stars in relation to 
the PMS tracks, that some of the distances chosen are not 
valid i.e. the star may not be associated with the chosen re¬ 
gion. Also, the spectrally determined Teff is more likely to 
be correct as it comes form spectra which has been directly 
observed from the star itself, opposed to a distance inferred 
from a possible association. A solution to this problem is cal¬ 
culating new distances to these outlier targets; ones which 
provide more sensible radii and agree with the spectrally 


determined temperature. To do this the stars are placed on 
the ZAMS at a point appropriate for their derived temper¬ 
ature; essentially, this is a lower limit to the luminosity of 
the star. This provides values of L*, i?*, M*, and an age. 
With the new ZAMS radii, revised distances are calculated 
from D/J?*. All objects affected by these ZAMS changes are 
noted in Tabled At this point all basic stellar parameters, 
relevant to this work, have been determined. 


4 BALMER EXCESS MEASUREMENTS 


With knowledge of the stellar parameters obtained for all 
targets a measurement of the Balmer Excess, ADb, can now 
be made. ADb is defined as the excess in flux above the 
intrinsic photospheric flux, seen across the Balmer Jump 
region (this region spans the wavelength range where the 
hydrogen Balmer series reaches its recombination limit ~ 
3640-3680 A). The UV-excess is weaker, in terms of energy, 
in lower mass stars, but is more readily visible due to their 
cooler photospheres, on top of which the excess can be seen. 
This UV-excess has been measured in both brown dwarfs 
jHerezeg_&^^iUenbrmi^ 20081 : iHerczeg. CnM_j£jIiljenbrandl 
20091; iRigliaco et al.l l2012l) and CTT s llCalvet et al.l liooil: 
Ordering et al. 2000l : Cal vet et al.l 12004 : llnglebv et alJ 
2 OI 3 II . From these past studies the current consensus to 
the origin of the excess is magnetospheric accretion. It 
has also been shown, in small samples, that an observable 
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Wovelength (A) 


Figure 3. Here is an example of step 2 in the stellar parame¬ 
ter determinations (see Section [3.2ll . where a reddened KC-model 
(black) is fitted to the observed photometry (red diamonds). 
The opposite is also shown, of dereddened photometry (blue plus 
signs) fitted to a KC-model with no reddening applied (dashed 
line). The level of reddening, Ay, is displayed in the figure. The 
fit provides a ratio between the distance to the star and its radius, 
D/i?*, as this is required to scale the model to fit the photom¬ 
etry. Also visible in this plot is how the C/-band magnitude is 
higher than the KC-model used, a possible indication of Balmer 
Excess. A clear IR excess can also been seen, starting at around 
the J-band; a typical feature of PMS stars. 



3.8 4.0 4.2 4.4 

Literature log(Temperature) [k] 


Figure 4. Shows the temperatures derived in this work in com¬ 
parison to literature estimates. The solid black line is the expected 
line of correlation and the dashed lines are a 0.1 dex deviation 
from this. The standard deviation between the two is only 0.02, 
with a mean offset of 0.03 towards this work. The dashed lines 
therefore encompass 3cr, showing the two samples are well corre¬ 
lated. The literature temperatures used and their references are 
provided in Appendix m 



3600 3800 4000 4200 4400 4600 

Wovelength (A) 

Figure 5. An example of the method 1 technique of measuring 
ADb is shown here. Initially the observed spectrum (dashed- 
black) and the intrinsic spectrum (red, a KC-model matching 
the spectral type of the target) have been normalised to 4000 A. 
Continuum points are selected from these between 4000—4600 A 
(shown as blue points). A ratio of these are provided in the bottom 
panel and they are fitted by a reddening law, which is extrapo¬ 
lated to 3600 A. This level of reddening correction is then applied 
to the original spectrum, with the result plotted in the top panel 
(grey). The SED of the model and the corrected spectra are now 
exactly the same between 4000-4600 A, allowing measurement of 
the Balmer Excess to be performed using only the region around 
3600 A. 


ADb in HAeBes stars can be explained within the same 
context (iMuzerolle et all 2004 Donel^w&g3rittainl I2OIII : 


iMendigutfa et al.l l2011bl : iPoeodin et al.ll2012l 'l . We aim to 
further our understanding of accretion in HAeBes by test¬ 
ing accretion within the context of MA to a large sample of 
HAeBes; this includes numerous HBe stars for which little 
investigation has been done. The Balmer Excess is defined 


ADb = {U- B)o -{U- H)dered (1) 

where, {U — B)o is the intrinsic colour of the target and 
{U — is the dereddened observed colour index. De¬ 

tailed below are the two best methods of measurement. 


4.1 Method 1 — Spectral Matching: Single Point 
Measurement 

The first approach to measuring ADb uses the spectral re¬ 
gion of the the UVB arm from 3500-4600 A, and adopts th e 
same techniques employed by iDonehew fc BrittainI ll201lli . 
This method requires the spectrum of the target to be com¬ 
pared against the intrinsic spectrum of a star of the same 
spectral type. The KC-models mentioned earlier are used 
here as the intrinsic star spectra. Following the calibration 
in Section [21 the spectrum of each target shows the correct, 
reddened spectral shape. This allows both the target and 
model spectra to be normalised to 4000 A, while preserving 
their spectral shape. Next, a correction for reddening present 
in the observed spectra is performed. To do this, the differ¬ 
ence between the measured continuum of the target and the 
model between 4000- 4600 A is fitted by a reddening l aw 
(the reddening law of iCardelli. Clayton fc MathisI (Il989l i is 
used here). This also provides a best-fit Ay - Extinction cor- 
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rection is applied to the whole spectrum, while maintaining 
the 4000 A pivot point for this correction. The resnlt of 
this method is that the spectral shape of the target is ad¬ 
justed such that the slope between the intrinsic model and 
the target spectrum now match. The success of this normal¬ 
isation is independent of the amount of extinction toward s 
the star (iMuzerolle et al.ll2004l : iDonehew fc BrittairJl2ojT]l . 
Fig [5] shows the application of this spectral slope matching 
technique along with an example output. 

To perform the measurement of ADs attention must be 
drawn back to Equation [T] where the magnitudes are now 
converted into a flux: 

( rpphot \ / jpdered \ 

7|h;Tj+2.51og(^^j (2) 


where F is the flux, with subscripts denoting the correspond¬ 
ing wavelength region, and the superscripts are: the intrin¬ 
sic flux - denoted ‘phot’, and the dereddened flux - denoted 
‘dered’. For these measurements the fluxes are monochro¬ 
matic. Now, consider the fact that the observed, dered¬ 
dened flux includes an accretion contribution, such that 
pAereA _ ^phot p^cc ^ This allows the above equation to 
be written as: 


ADb = 2.5 log 


F phot 

_u _ 

T-iphot 
^ B 


I TTiacc 

+ i'u 

-l-^acc 



(3) 


This equation can be reduced through the use of a normali¬ 
sation factor a"”™, where -f x a"”™ = Fl^°\ 

This normalisation across the B-band is performed automat¬ 
ically by matching the slope of the spectrum of the target 
to the intrinsic spectrum’s slope (see the steps mentioned 
earlier). In essence 0 “°””“ represents a reddening law. This 
gives us the final form of the ADb equation: 


ADs = 2.5 log 


F phot,norm 

u 


+ F, 


acc,norm 


-iphot 


(4) 


By these definitions, the _|_ ^acc.norm 

flux observed from the target spectra, and can be 

taken from a KC-model of the same spectral type. Since 
the spectrum obtained is of medium resolution we adopt a 
narrow, monochromatic, range over a typical broadband fil¬ 
ter to represent the 17-band magnitude. This also gives us 
better precision in measurements. The wavelength region of 
measurement is 3500-3680 A. This is chosen as it is beyond 
the Balmer recombination limit. However, two of the echelle 
orders of X-Shooter overlap in this region, and the SNR in 
an echelle order decreases as wavelength decreases. There¬ 
fore, to minimise errors, the 3500-3600 A region from echelle 
order 21 and the 3600-3680 A region from echelle order 20 
are measured and combined to give the most accurate result. 


4.2 Method 2 — B-Band Normalised, Multi-point 
Measurements 

An alternate method of measuring ADb is given by 
iMendigutfa et ahl (l2013l l. which also does not require the 
reddening towards a star to be known. This method covers 
a larger wavelength range, requiring measurements of both 
the 17-band and F-band points. These two points are mea¬ 
sured from the observed spectra and a KC-model of the same 


spectral type (the same model as in method 1), after nor¬ 
malisation to the B-band. Rather than correcting for Ay, as 
in the previous method, reddening independence is achieved 
by expanding Equation [T] and substituting in an expression 
for each reddening component: Ax = Avikx/kv), where Ax 
and Av are the extinction at any given wavelength and in 
the F-band, respectively. Similarly, kx and ky are the opac¬ 
ities for any given wavelength and the F-band, respectively. 
Applying the expression for Ay to Equation [T] gives: 

ADb = {U- B)“ - (17 - + Ay(!p^- (5) 

\ky ky J 

The superscript ‘inf refers to the intrinsic magnitudes (from 
a KC-model in this case) while the superscript ‘obs’ refers to 
the observed magnitudes (from the observed spectra). The 
values of the opacities are deter mined by the reddening law 
adopt ed. The reddening law of ICardelli, Clayton fc MathisI 
(Il9891) is used here with an Rv=3.1, providing ku/ky = 1.57 
and ks/ky = 1.33. To remove the Ay term the rela¬ 
tionship between Ay and colour excess needs to be used: 
Ay = IlyE{B — F). At this point it should be noted that 
the method is now reddening independent, since Ay has 
been removed, but remains dependent on the reddening law 
adopted, as this affects the opacity ratios. This new form for 
the Balmer Excess is: 

ADb = {U - - (17 - B)°'’=-f 

[{B - F)°‘’= - (B - F)“] (6) 

\ky ky J 

which can be expressed in terms of flux, instead of magni¬ 
tudes, as follows: 


ADb = 2.5 log 


2.5 Rv 


F obsnorm mphot 
1 mphot 


I 


-h 


kif kB 


k\ 


log 

7 V F^^a'^ 


^phot 


^phot 


(7) 


where 0 :"°'^™ has been added and is a normalising factor for 
the B-band, as seen in method 1, but the normalisation is 
instead performed such that the spectra will be unity at 
4400 A. Note that all these fluxes are considered monochro¬ 
matic, with centres at the usual Johnson UBV wavelengths. 
This normalisation allows the equation to reduce to its final 
form: 


ADb = 2.5 log 


^phot 


2.5Rv('^-^ 

^ kv ky 


log 


nobs,norm 

_V _ 

^phot 


(8) 


in this form it can be seen that only four points need to 
be measured in order to obtain ADb (two from the target 
spectra, two from the model). 


4.3 Comparisons and Checks 

The two methods used are similar but have some subtle dif¬ 
ferences. One, is that the central wavelength for the B-band 
normalisation is different between the two; it is centred at 
4000 A for method 1, and is centred at 4400 A for method 2. 
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Ay from Photometry Fitting 


Figure 6 . This shows a comparison of a pseudo-Fly, extracted 
from the method f measurement of ADg depending upon how 
much the spectra was adjusted, against the Ay determined from 
the photometry fitting in Section 13.21 The solid black line is the 
line of correlation, while the dashed lines are fcr deviations, of 
0.60 mag from this. Only 8% of the targets are outside 2 (t; these 
ones often have the largest Ay values (these will be discussed in 
the text). Errors in the photometric Ay are typically 0.05—0.15 
(about the width of the points), and are given in Table [J] 


The next difference is that method 1 performs a reddening 
correction using a section of the observed spectrum and re¬ 
lies on matching this to a stellar model. On the other-hand, 
method 2 avoids having to make a reddening correction by 
incorporating the adopted reddening law into the equation 
for ADb, and applies this over a much larger spectral region. 
Also, both approaches have been adapted from a definition 
which was based on broadband photometry. Therefore, some 
checks need to be made to see whether both approaches are 
comparable to each other. 

The first check is between how the Ay values deter¬ 
mined in Section O compare with the Ay values extracted 
from method 1; as the fitting between 4000-4600 A can be 
used to infer an Ay value. Figure[6]displays this comparison. 
In the figure the standard deviation between the two is found 
to be 0.60 mag, and is represented by the dashed black lines. 
Within this la interval 79% of the sample are included. This 
helps to highlight that the majority of the sample are tightly 
correlated while the outliers are more extreme and actually 
skew the standard deviation towards them. There are 7 stars 
showing differences greater than 2a from the mean. One of 
these is VY Mon which has the lowest SNR of the objects in 
the blue because it is very extinct. This makes the spectral 
shape adjustment more difficult and less accurate than other 
targets. The other outliers often have large Ay and/or large 
ADb values. This is not entirely unexpected as a signifi¬ 
cant excess can affect the SED shape of the spectra, which 
would complicate both photometry fitting and the spectral 



ADg Method 2 

Figure 7. A comparison is made between the two different meth¬ 
ods of measuring ADg (detailed in Section |4.II A line of 
expected 1:1 correlation is shown in black. For clarity individual 
error bars are not plotted due to the tightness of the points. In¬ 
stead, a typical error bar of 0.04 mag is plotted in the bottom 
right corner. The actual errors range between 0.04—0.10 mag, and 
are provided in Table [3] 


shape adjustments performed. In general, for HAeBes this 
is less likely as they are already very hot and the excesses 
need to be very strong to significantly affect the SED. One 
source of discrepancy lies in how the photometric method 
is coarse but covers the BVRI points, while the spectral 
method covers a very narrow wavelength range of 4000 - 
4600 A, but with a greater accuracy in that region. The 
photometry used is also not simultaneous with the spectra; 
variability could therefore also play a role in the differences. 
Ultimately, this scatter is quite low with few outliers; this is 
more than acceptable considering the above factors and the 
standard reddening law adopted in both cases. 

The next check is to see how ADb varies between the 
two methods of measurement; Figure [7] shows the compari¬ 
son. There is a systematic offset of ~ 0.02 towards method 
2 producing higher values, while the standard deviation of 
scatter between the two methods is 0.04 mag. These dif¬ 
ferences are less than the systematic error on measurements. 
Since the original ADb equation, Equation[T] can be seen to 
contain a dereddened term, the differences can be mostly at¬ 
tributed to how the reddening corrections are made in each 
case; though the normalisation of the spectra and points of 
measurement also influence the result. The method 1 covers 
a small wavelength range of 3600 - 4600 A, of which only the 
4000 - 4600 A region is used for the reddening correction. 
This means that this approach is not particularly sensitive 
to a given reddening law due to the small wavelength range 
it covers, and can be deemed reddening independent for low 
levels of extinction {Ay < 10). On the other hand, method 2 
depends more upon the adopted reddening law than method 
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ADg from Previous Works 


Figure 8. A comparison is drawn here between the final ADb 
measured in this w ork versus the ADb measur e d by other authors 
in the literature llDonehew &: BrittainI I2OIII : iMendigutia et ^ 
bnild: IPogodin et a 1J I2012I1 . The difference in temperature be¬ 
tween the two sources is calculated as a percentage of the total 
stellar temperature (cooler temperatures than the literature are 
white, while hotter ones are black). The size of each symbol re¬ 
flects the difference in log(g) measured. Overall the largest devi¬ 
ations in ADg are for the objects with the greatest differences in 
stellar parameters. 


1, because it covers a larger wavelength region of 3600 - 
5500 A. Depending on the Rv selected the resulting opacity 
ratios, seen in Equation can change substantially, which 
in turn alters the measured ADb- Changing Rv in method 
1 does not noticeably affect ADb, for low Av values, as 
it is always the spectral profiles which are being matched. 
Through this matching the Av used will change to retain 
the SED shape and keeps ADb the same. Returning to the 
figure, a few outliers can be seen between the two methods; 
the majority of these are objects with high extinction, or 
which were identified as having a discrepant Ay between 
the photometric method and the spectral method in which 
they were determined. 

Overall, consistency is apparent between the methods 
employed here as the majority of measurements from each 
method lie within the errors of each other (see Table O). 
Based on the above analysis, we deem the methods equiv¬ 
alent. Therefore, in each case an average of the two will be 
taken for the final result; unless one method has a lower mea¬ 
surement error, then that method will be favoured over the 
other (this can occur depending on emission lines in both the 
measurement and normalisation regions). The ADb value 
for each star, along with the errors and method(s) used to 
obtain it, are detailed in Table [S] The errors given in ADb 
appear large when compared with the value of ADb itself. 
It should be noted that the detections are above Su and the 
enhanced errors are mostly due to taking the logarithm of 


a ratio, see Equation [T] where an error of 1% in the con¬ 
tinuum detection can translate to more than a 30% error in 
ADb (depending upon how small the difference is between 
the intrinsic and observed spectra). 

A comparison of the ADb values determined in this 
work versus previous values published in the literature is 
shown in Figure (S] as a consistency check. The majority of 
the measurements are clustered at values < 0.4 mag, with 
literature values showing a slightly larger spread in ADb 
than our sample. The main source of deviation between this 
work and the literature can be attributed to the Teff and 
log(g) parameters used for each star; as these differ so will 
the intrinsic spectra form which ADb is measured. The fig¬ 
ure shows this clearly with a number of objects having de¬ 
viant ADb and stellar parameters, where the largest varia¬ 
tions in ADb are indeed the stars with the largest changes 
in Tefi and log(g), compared to the literature values. How¬ 
ever, there is one star whose deviation in ADb cannot be 
explained by the changes in Tefi and log(g) alone. Instead, 
the deviations may also be compounded by genuine variabil¬ 
ity of the star and/or accretion rate. Such variability can be 
seen within the lite r ature, and in single st ars themselves 
dPogodin et ahllioil : iMendieutfa et al.l]201^ 1. Additionally, 
intrinsic features in the spectra can contribute to differences 
too; the approach in this manuscript uses monochromatic 
points from spectra, whereas the majority of comparison 
stars primarily use broad-band photometry. Overall, the ma¬ 
jority of sources are in common, within the errors, and most 
discrepancies can be explained by the adoption of stellar 
parameters. 


5 ACCRETION RATES 

Accretion rates are an important parameter of pre-main se¬ 
quence stars. They provide an insight into how the stars are 
evolving, along with the impact this will have on disc-star 
interactions, and may even have repercussions on planet for¬ 
mation. 

5.1 Magnetospheric Modelling 

In this work the measured ADb is used to calculate Macc 
using accretion shock-modelling within the context of MA. 
This theory is adopted in order to test its applicability 
to a wide sample of HAeBes. The main assumption here 
is that the excess flux visible over the Balmer Jump re¬ 
gion is produced by shocked emission from an in-falling 
accretion column. A detailed description of the magneto- 
spherical ly driven accretion column and shock-modelling is 
given bv ICalvet fc Gullbrind ll 19981 . hereafter CG98), while 
a description of its ap plica tion to HAeBe stars is gi ven in 
iMuzerolIe et al.l (12004 ) and lMendieutfa et al.l (l2011bi) . Here 
we summarise the key points of those papers and detail how 
they work in regard to this sample. 

Firstly, the magnetic field lines of the star interact with 
the disc and truncate it at the truncation radius, Ri. It is 
generally accepted that the truncat ion radius is close to, or 
insid e, the co-rotation radius, i?cor, (lKoenig]||l99ll : IShu et al.l 
Il994l . CG98). For this work Ri is chosen to be 2.5 7?*, as this 
has been s hown to be an appropriate value which is often les s 
than Rcor llMuzerolle et al.l[20041 : [Mendieutfa et aLllioilbl) . 
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Table 3. Table of Accretion Rates. Column 1 gives the target name. Columns 2-4 give the measured 
Balmer Excess, filling factor and derived accretion rate. ADb errors are rounded to the closest 0.01 
and include all systematic errors too. Column 5 details by which method the values were obtained. 
Column 6 gives the accretion luminosity. Finally, Column 7 notes which stars can have their excess 
modelled successfully by MA. 


Name 

ADb 

/ 

log(Afacc) 

Method(s) 


log(Lacc ) 

Achievable 


(mag) 

(%) 

[M 0 /yr] 

Used 


[L©] 

by MA 

UX Ori 

^ 0.04 

^ 0.7 

^ -7.26 

Method 1 & 

2 

^ 0.13 

y 

PDS 174 

^ 0.02 

s£ 2.8 

^ -6.76 

Method 2 


^ 0.92 

y 
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u.iy_Q Q 5 


* ■^'^-0.28 

Method 1 & 

2 

u.oo_Q 32 

y 

HD 34282 

05 

1 7 + 0.9 
-*^•'- 0.8 

-7 69"*"*^'^^ 

+ 'jy_0.59 

Method 1 


y-yO-0.61 

y 

HD 287823 

Q 1 K-h0.05 
U- 10 _o 05 

O.U_o 7 

7 1 q + 0.18 

‘ •-^'^-0.23 

Method 1 & 

2 

0 o-t+ 0.20 

'^•'^'-0.24 

y 

HD 287841 

C 0.05 

^ 0.8 

^ -7.82 

Method 1 & 

2 

^ -0.32 

y 

HD 290409 

C 0.07 

^ 2.1 

^ -7.31 

Method 1 & 

2 

C 0.25 

y 

HD 35929 

y-J-y_o.05 

1 0 +°'^ 

-6 37+*^'i^ 
-0.26 

Method 1 & 

2 

0.87l°;®° 

y 

HD 290500 

'^•^-'--0.05 

’^•-‘^- 1.5 

'^•-'--‘--0.17 

Method 1 & 

2 

l.29l°;®® 

y 

HD 244314 

f) 12 + 00 ^ 

U-J-^_0.05 

2 4+0-’^ 

7 1 9 + 0.20 
■'•1^-0.25 

Method 1 & 

2 

U.OO_Q 26 

y 

HK Ori 

U-00_o.o5 

97 7 + 6.0 

6 1 7+0-17 
-O-i<- 0 .i 6 

Method 1 & 

2 

1 qq+0.19 

1 . 00 _Q ig 

y 

HD 244604 

U-UO_0.05 


7 Qq+0.26 

-i .zz_Q 32 

Method 1 


r\ qQ + 0.28 
U.zz_o 34 

y 

UY Ori 

C 0.02 

^ 0.6 

^ -7.92 

Method 1 & 

2 

^ -0.35 

y 

HD 245185 

C 0.07 

^ 2.3 

^ -7.29 

Method 1 & 

2 

C 0.29 

y 

T Ori 

C 0.05 

^ 1.0 

^ -6.54 

Method 1 & 

2 

C 0.79 

y 

V380 Ori 

n q’7”I"0-05 

O.o^_o.o5 

80.311®;" 

-5.34l°;l® 

Method 1 & 

2 

9 1 Q+0 .31 
"^•-*-"^-0.16 

y 

HD 37258 

U.i^_0.05 

4-51M 

-6.9810-14 

Method 1 & 

2 

0.58l°;4 

y 

HD 290770 

Q 1 k-FO.05 
U- 10 _o 05 


-6 74+*^-i^ 
'^-0.14 

Method 1 & 

2 

0.82l°;4 

y 

BF Ori 

Q 1 k-FO. 05 
U-J-O_o.o5 

3-6lg-| 

D.DO_o 25 

Method 2 


0.77l°;4 

y 

HD 37357 

U-oU_o.o5 

lo.+M 

6 42+*^-'^® 

-O-+^_0.06 

Method 1 & 

2 

l-08lg:g^ 

y 

HD 290764 

0 21 

o C+O.S 
o.o_o 7 

D.OD_q 22 

Method 1 & 

2 

0-80l°;®i 

y 

HD 37411 

Q 1 k-FO.05 
U- 10 _o 05 

4 q+1'3 

7 1 q+0.24 

< •-LO_o 34 

Method 1 & 

2 

0 47+0-29 
_0.38 

y 

V599 Ori 

^ 0.01 

^ 0.1 

^ -7.67 

Method 2 


C -0.33 

y 

V350 Ori 

n 1 c-t-0.05 
U- 10 _o 05 

q 7 + 1.0 

'^•^-0.9 

-6.95l0;4 

Method 1 & 

2 

0 >55+0-10 

U.OO_Q 24 

y 

HD 250550 

y-'jy_0.05 

-*-'•-*--3.1 

-5.63l°;4 

Method 1 & 

2 

1 09 + 0.32 

l-o2-o 17 

y 

V791 Mon 

y-j-y_o.o5 

97 k + 10.2 

-6 16+*^'2^ 
o-io_o.35 

Method 1 & 

2 

1 ck+0.32 
1-3^-0.39 

y 

PDS 124 

nil +0.05 
^•-*--*^-0.05 

4 1 +1-3 

71-1 +0.13 
- + 11-0.19 

Method 2 


U. 0 U_Q 20 

y 

LkHa 339 

f) 1 O+0.05 
U-J-O_o.o5 

K 0 + 1.5 

0 . 0 _ 2 .i 

-6.8ilg:4 

Method 2 


0 75+0-13 
^^+^- 0.22 

y 

VY Mon 


1 ‘ -^-13.0 

-5.50l0;g 

Method 1 & 

2 

i-yo_o .66 

y 

R Mon 

0-861^5 

- 

- 

Method 1 & 

2 

- 

n 

V590 Mon 

- 

- 

- 

- 


- 

- 

PDS 24 

^ 0.05 

^ 1.9 

^ -7.25 

Method 2 


^ 0.31 

y 

PDS 130 

y--LO_o.o5 

6-6ll®5 


Method 1 & 

2 

1 qq-I- 0.17 
i-^^-0.l5 

y 

PDS 229N 

y-yy_o.o5 


-6-67l0;4 

Method 2 


y-yO-o 36 

y 

GU CMa 

0 14+005 
U.14_o 05 

60.3l?|'l 

-5.0010;4 

Method 1 & 

2 

9 79+0.41 

^+^-0.36 

y 

HT CMa 

nil +0.05 
^•-‘--‘^-0.05 

4.311® 

6 61 +*^- 1 ® 
-D.D 1 _q 19 

Method 1 & 

2 

0.89l®;4 

y 

Z CMa 

1-O5_o.o5 

48.011"/ 

-o.Ui_Q 19 

Method 1 & 

2 

4.05l®;i 

y 

HU CMa 

0 14+005 

U.l+_o.o5 


0 oc+ 0.11 

O.OO-O 22 

Method 1 & 

2 

1 97 + 0.11 
-^•^'- 0.22 

y 

HD 53367 

y--Ly_o.05 

- 

- 

Method 1 & 

2 

- 

n 

PDS 241 

0-051^5 

21 . 61 ^ 3 ® 

^•3O_0.06 

Method 1 & 

2 

9 9 ^+ 0 .37 

z.^o_Q 23 

y 

NX Pup 

0.081^5 

0 Q+0-4 

-6.96l0;2i 

Method 2 


0 34+0-23 

y-o^_0.34 

y 

PDS 27 

0 1 7+013 

40 0+35.0 
4U.U_3g 0 

-3.9610-JS 

Method 1 & 

2 

3 40+0-78 

o.+y_o.82 

y 

PDS 133 

l-^0_o.o5 

- 

- 

Method 1 & 

2 

- 

n 

HD 59319 

U-‘JO_0.05 

0 4+2.3 

-5.76l0;4 

Method 1 & 

2 

1 65+0-26 
1-OO-0.15 

y 
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Table 3 - continued 


Name 

ADb 

/ 

log(Macc) 

Method (s) 


log(i/acc) 

Achievable 


(mag) 

(%) 

[Mo/yr] 

Used 


[L©] 

by MA 

PDS 134 

^ 0.03 

s: 3.0 

s£ -5.60 

Method 2 


^ 1.82 

y 

HD 68695 

'-'•U0_o 05 

1 q+0.8 

-‘^•'^-0.1 

Y '7Q+0.30 
- ^ /5_o_38 

Method 2 


-0 17+6-34 
1^-4'-0.41 

y 

HD 72106 

Q O1+0.05 
‘J-ol_o.o5 

7-+1:3 

-6.2i1°;13 

Method 1 & 

2 

i.20l°;ig 

y 

TYC 8581-2002-1 

Q 1 r+0.05 

‘J--LO_0.05 


-6.58l°;i° 

Method 1 h. 

2 

0.88l°;“ 

y 

PDS 33 

^ 0.04 

1.2 

^ -7.84 

Method 1 h. 

2 

^ -0.21 

y 

HD 76534 

^ 0.01 

s; 1.7 

^ -6.95 

Method 1 h. 

2 

^ 0.77 

y 

PDS 281 

- 

- 

- 

- 


- 

- 

PDS 286 

U-U/_o.o5 

64 fi+35.4 
0^-0_39,6 

5 41 +0-20 
-0.4i_Q gg 

Method 1 h. 

2 

9 55+0.62 

z.oo_Q 

y 

PDS 297 

^ 0.01 

s: 0.4 

^ -7.60 

Method 2 


^ -0.12 

y 

HD 85567 

^•oo_0 05 

- 

- 

Method 1 h. 

2 

- 

n 

HD 87403 

0 05+0 05 
‘J-U0_o 05 

1 p:+0.9 
4-3_o.2 

5 82+*^'^*^ 
^•o^_0.07 

Method 2 


1 48+*^'^^ 

y 

PDS 37 

0 10+0-22 

4U.U_3g Q 

3 50+0-60 
o.oo_i g2 

Method 1 


q oc;+0.65 

0.00_1 g2 

y 

HD 305298 

'J-^O_0.05 

- 

- 

Method 2 


- 

n 

HD 94509 

- 

- 

- 

- 


- 

- 

HD 95881 

^ 0.05 

s: 1.5 

^ -5.65 

Method 1 & 

2 

^ 1.63 

y 

HD 96042 

0 12+005 

Q4 4+5.6 
y+'+_38.5 

4 57+0-03 
^•^'-0.28 

Method 1 & 

2 

q 1 q + 0.32 
o.iO_0 30 

y 

HD 97048 

^ 0.01 

sC 0.4 

-8.16 

Method 2 


^ -0.55 

y 

HD 98922 

^ 0.01 

s; 0.4 

^ -6.97 

Method 2 


^ 0.41 

y 

HD 100453 

^ 0.01 

s; 0.1 

^ -8.31 

Method 1 & 

2 

^ -0.92 

y 

HD 100546 

n 1 q+O- 05 
0.1o_0 05 

6 1+1® 

7 04+*^'^^ 

-^.U4_q 

Method 1 & 

2 

0.56l°-14 

y 

HD 101412 

^ 0.04 

s: 1.2 

^ -7.61 

Method 1 & 

2 

^ -0.04 

y 

PDS 344 

^ 0.03 

s: 3.5 

^ -7.02 

Method 1 & 

2 

^ 0.68 

y 

HD 104237 

0 17+005 

9 q+0.7 
^■*- 0.6 

-6.68l°;i3 

Method 1 h. 

2 

0 70+6 -16 

6 ’- ' 6 '_o .22 

y 

V1028 Cen 

0 10+005 

'J--L»J_o.05 

10.613-7 

-5.76l°;i® 

Method 1 


1 76+6-^6 
-*-• '^-0.24 

y 

PDS 361S 

0 12+005 
0'--L2_o 05 

26.21®;® 

5 26+^’^^ 
^•^0'-0.20 

Method 1 h. 

2 

9 qc;+0.27 
^.o<j_o 23 

y 

HD 114981 

^ 0.06 

s: 8.1 

^ -5.48 

Method 1 & 

2 

^ 2.03 

y 

PDS 364 

n 9Q+0.05 
0.^o_0.05 

26.81®;^ 

-6.05l°;i3 

Method 1 h. 

2 

1 50+0.13 
1.0O_Qj^2 

y 

PDS 69 

0 O1+0.05 

0).ol_o.o5 

62.81®®;® 

5 Q9+0.21 

-O.OZ_Q 21 

Method 1 h. 

2 

9 Qq + 0.30 
^•-iO_0.27 

y 

DG Cir 

0 7Q+0 05 

o. ^y_o.05 

- 

- 

Method 1 h. 

2 

- 

n 

HD 132947 

0 06+0 05 
0'-UO_o 05 


-6 71+0-17 
' ^-0.42 

Method 1 & 

2 

0 73+0-18 
''^-0.42 

y 

HD 135344B 

0 07+0 05 
Oi-U/ _o 05 

0 7+0-3 

O'- '-0.3 

<J q-7'+0.24 
■'•'^'-0.41 

Method 1 & 

2 

0'-O'i_o.42 

y 

HD 139614 

o'-uy_o 05 

1 k+ 0.6 

-^•^-0.5 

-7 63+0-20 

{ .UO_Q 3Q 

Method 1 


0'-i0»_0.31 

y 

PDS 144S 

^ 0.01 

C 0.1 

^ -8.35 

Method 1 


^ -0.90 

y 

HD 141569 

n 05+0 05 
0'-U0_o 05 

1 5+0.9 

7 65+0-33 

+03-0.47 

Method 1 


0 05+*^'^^ 

U- 0 iO_Q 50 

y 

HD 141926 

0 20+0 05 
OI-^U_o 05 

- 

- 

Method 1 & 

2 

- 

n 

HD 142666 

^ 0.01 

^ 0.1 

^ -8.38 

Method 1 & 

2 

^ -0.93 

y 

HD 142527 

0 06+0 05 
‘J.UD_o 05 

0 0+0-3 

O'-'^-O.S 

7 45+0-19 

+ 43_q 4 g 

Method 1 


0 09+*^'^® 
u.uy_o.48 

y 

HD 144432 

0 07+0 05 
Oi-U/ _o 05 

1 0+0-5 

^7 qo + 0.22 
- i •30_q 4 Q 

Method 1 & 

2 

0 02 +*^'^^ 
6'-6'^-0.41 

y 

HD 144668 

0 '-^oi_o 05 

o q+0.9 
'^•^-0.8 

6 25+^+® 

o-2O_0 19 

Method 1 & 

2 

110 + 6+6 
-*- + 6 '- 0.22 

y 

HD 145718 

^ 0.01 

C 0.2 

^ -8.51 

Method 1 


^ - 1.01 

y 

PDS 415N 

^ 0.04 

^ 0.5 

^ -8.45 

Method 1 h. 

2 

^ -0.91 

y 

HD 150193 

0 07+0 05 
Oi-U/ _o 05 

-‘-•'^-0.7 

-7 45+'^'^5 

+^3_o 43 

Method 2 


0 10 + 6-^6 
0'-10_o.44 

y 

AK Sco 

^ 0.04 

^ 0.4 

^ -7.90 

Method 1 


^ -0.52 

y 

PDS 431 

0 11+0.05 

4 3+^'^ 

O.UO_Q 22 

Method 1 h. 

2 

1 34+6-21 
-*-•0^-0.24 

y 

KK Oph 

^ 0.05 

^ 1.0 

^ -7.84 

Method 2 


^ -0.29 

y 

HD 163296 

0 07+0 05 
‘J-U/_o 05 

1 O+ 0.8 

-*-•^- 0.8 

7 40+0-14 

+4y_o.3o 

Method 2 


0.08l°;“ 

y 

MWC 297 

oil +0.08 

‘J.il_o.08 

56 ‘^+‘^5.7 
ou.o_26.5 

5 16+0-25 
-0.1D_q 43 

Method 1 h. 

2 

2-62lg;tg 

y 
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Figure 10. This figure shows KC-model atmospheres (red) for a 10000 K star, on the left, and a 20000 K star, on the right. In both 
cases .^=10^^ erg cm“^ and log(Macc)=-6.5. These allow the accretion flux (black, dashed) to be calculated. Adding this flux to 

the intrinsic photosphere provides the observed spectra (black, dotted). This spectra is then normalised to 4000 A (grey), so that its 
continuum matches the intrinsic between 4000—4600 A, allowing ADg to be measured via the method 1 approach. The resulting ADg 
values are given in the plot, demonstrating how they vary depending on the temperature of the star. 


Rcoi can be smaller than the adopted 2.5 i?*, as is the case is 
for fast rotators, but this will not affect the derived accretion 
rate significantly i.e. for a very small truncation radius of 
Ri = 1.5 R© the resulting accretion rate would be less than 
a factor of two different from one where Ri =2.5 R©. 

At the truncation radius material is funnelled by the 
field lines and falls at speeds close to free-fall towards the 
stellar surface, where it shocks the photosphere upon im¬ 
pact. The velocity of the infalling material, Va, is given as: 






(9) 


The velocity can be related to the accretion rate via the 
density. This is because Macc is flowing at the same rate as 
the velocity through an accretion column, which also cov¬ 
ers a given area of the star. Therefore, the density can be 
expressed as: 


P = 


M,, 


AVa 


( 10 ) 


where A is the area of the star covered by the accretion 
column, defined as A = /47rRj, and / is a filling factor such 
that / = 0.1 would be 10% surface coverage. The filling 
factor is required as we consider the accretion to be funnelled 
though a column, rather than being evenly distributed over 
the entire stellar surface. Putting this in terms of energy, 
the total inward flux of energy of the accretion column is: 

= {l/2)pV^a ( 11 ) 


This amount of energy is carried into t he colu mn and must 
be re-emitted back out of the star (see ICGOSi for details on 
this energy balance). This means the total luminosity from 


the accretion column, as given in 




can be written as: 


Lcoi = {^ + F.)A = C + f.A = CLacc + F.A 

( 12 ) 

where F^, is the intrinsic flux of the stellar photosphere. Lace 
is the accretion luminosity, and ^ = 1 — (R*/i?i)- The accre¬ 
tion luminosity is defined as Lacc = G'MaccM*/i?*. 

As shown in iMendigutfa et al.l ll2011bl ). the column lu¬ 
minosity is Lcoi = LcoiA, where Led is the flux produced 
by the accretion column. This total amount of flux can be 
expressed as a blackbody function, where Lcoi = Sim¬ 

ilarly the same can be done for the photosphere, L* = aTf 
This results in oTf. 

At this point the unknowns are /, Tcoi, and Ma.cc- 
Tcoi has just been shown to be governed by the amount of 
energy flowing onto the photosphere, and by the tem¬ 
perature of the photosphere itself, T*. For each star Tcoi is 
determined using the temperatures we derived, and by fix¬ 
ing — 10^^ erg cm“^ as this has been shown to 

provide appropriate filling factors of ^ 0.15 in the major¬ 
ity of cases in HAeBes st udied so far (iMuzerolle et al.ll2004l : 
IMendigutfa et aLll2011bl) . This leaves only / and Mace re¬ 
maining. Ma.cc can be determined from ADg by making use 
of the equations above; for which there is a unique ADg vs. 
Macc combination for each star due to its stellar parame¬ 
ters. To obtain this curve, Macc values are tested between 
10“^-10“^°M©yr“^. With fixed, and all the other stellar 
parameters known, the filling factor corresponding to each 
Macc value is found through the following equation (which 
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Figure 9. The relationship between ADb and Ma.cc is shown 
here, for a series of stars with different temperatures (labelled in 
the legend). The stellar parameters used for each case are ones 
typical for a ZAMS star of the temperature in question. It is 
apparent that the same ADb value would result in a higher ac¬ 
cretion rate in hotter stars. The filled circles display the point 
at which / =0.1, where the accretion column covers 10% of the 
surface. Similarly, the filled diamonds are where / =1.0 (full cov¬ 
erage). In all cases ^=10^^ erg cm~^ A“^. 


is a rearrangement of the 2^^^ and 3^^^ terms in Eauation ll2p : 


/ = c 


/GMaccM* 

( R. 




(13) 


The Tcoi determined previously is used to make a black- 
body, which represents the accretion hotspot, and multiply¬ 
ing this by / gives the excess flux. The excess flux is then 
combined with a KC-model, determined using the relevant 
stellar parameters. From this, ADb can be measured. This 
is repeated for all Ma.cc and / combinations. The result pro¬ 
vides a unique ADb vs. Macc curve, which the accretion 
rate can be read from. 

Figure [5] gives the ADb vs. Macc curves for a series 
of different temperature stars (for simplicity in the figure 
their other parameters are taken from the ZAMS). The fig¬ 
ure demonstrates how the same ADb, measured in two dif¬ 
ferent temperature stars, can refer to wildly differing accre¬ 
tion rates. Also, the T=10000 K curve has the highest ADb 
value as the size of the Balmer Jump peaks at around this 
temperature. 

Figure fTUl demonstrates the same concept of ADb vs. 
Macc changing as a function of temperature, as shown by 
the curves in Figured Although, this figure also highlights 
how the excess flux impacts the appearance of the spectra 
too. There are two cases in the figure, one for a star of 
10000 K, and the other for a star of 20000 K. It can be seen 
for Macc = 10 “®'® that the resulting ADb changes from 0.41 
for the 10000 K star, to only 0.17 for the 20000 K star. This 
is why the calculation of separate ADb vs. Macc curves, for 


Figure 11. Macc versus T^ff is shown for each object, along with 
literature values for a comparison. In general it can be seen that 
Macc increases with temperature, with a scatter of 2-3 orders of 
magnitude in Macc throughout. 


each star, are crucial. It also demonstrates that the SED 
shape is not significantly affected by the excess longwards 
of 4000 A, which means that the approach of methods 1 
and 2 remain valid (as does the photometry fitting). Table 
[3] contains the Macc values for each star using an individual 
curve for each star. 


5.2 Literature Comparisons 

Comparisons of Macc derived in this work are made against 
previous detections in HAeBes and CTTs. In the first com¬ 
parison, Figure [11] places this sample against other stars 
from the literature in which Macc has also been determined 
directly using ADb - For the HAes, ^10000 K and lower, the 
range in this work appears similar to previous works, with 
Macc spanning from anywhere between 10 “®- 10 ~® MQyr“^, 
with the exception of one star at ~ 10~® M 0 yr“^ (Z CMa, 
which is likely a very young HBe star based on it’s mass 
of 11 Mq ). The HBes closest to the HAes show a simi¬ 
lar range in magnitude of 10~^-10“^ M 0 yr“^. The scatter 
then decreases once the temperature has increased beyond 
20000 K, where Macc spans lO'^-lO”"* M 0 yr“^. This de¬ 
crease can be partially attributed to a detection effect, as 
the temperature of the star increases the observable ADb 
will decrease. Therefore, if the temperature of the star is 
very high then low accretion rates will be undetectable via 
the Balmer Excess method. This is supported by the ADb 
vs. log(Macc) curves in Figure (Hj Returning to Figure [TT] 
comparisons are als o drawn against previous ly published 
accretion rates. The iMendieutfa et all ll2011bl l sample has 
a slightly larger scatter showing some Macc detections be¬ 
low our findings, this can again be attributed to detection 
limits in this work. But it can also be seen that there are 
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Figure 12. Macc versus M* is shown for each star, where possi¬ 
ble, alo ng with additional HAeBe sources from |Mendi g|utfa et ^ 
ll2011]Jl . and CTTs from iNhtta. Testi Randichl ||2006^ ■ A red 

solid line fit to all of the points, excluding upper limits, is shown 
J\/T ^ 71 T2.47zb0.07 
of Macc OC 



Figure 13. Lace versus L* is shown for each star in this work, 
along w ith additional HAeBe sources from fM endig utia et al] 
ll2011bl l. and CTTs from iNatta. Testi &: Randichl (120061 ') . A best 
fit is obtained of Lace oc which is plotted as a solid 

red line and excludes the upper limits. 


many stars in iMendieutfa et al.l ll2011bli which have accre¬ 
tion rates about an order of magnitude higher than our find¬ 
ings. The exact reason for the discrepancies is unknown, but 
it is likely to be a combination of the two different types of 
dataset, spectra and photometry, and the different methods 
of measurement used i.e. the photometric method requires 
dereddening to be performed prior to measurement of ADs. 
Variability may also play a role. 

Comparing ou r res ults with the work of 
iDonehew &: BrittainI (1201111 we find a systematically 
higher accretion rate for objects hotter than 10000 K, 
the HBes, of around 1-2 orders of magnitude. This can 
be attributed to their adoption of a single ADb vs. Macc 
relationship for all of their objects. Whereas in this work, 
the relationship between the two has been calculated 
on an individual basis for each star, based on its stellar 
parameters (see Figure [9ll. Therefore, they are not directly 
comparable. 

A comparison is made of Macc vs. M* in Figure 1121 
which includes literature values too. More specifically, this 
comparison looks at how the results of this sa mple compare 
to the HAeBes from lMendieutfa et akl ll2 011bf), along with a 
look a t lower luminosity CTTs from lNatta. Testi fc Randichl 
ll2006ll . A trend is seen of increasing accretion rate with 
increasing stellar mass; the fit shown in the figure gives 
Macc oc Af2.47±o.07 position of the HAeBes obtained 
in this work show agreem ent with the values obtained from 
IMendieutfa et akl (l2011bli . However, it can be seen at around 
the HAe mass range, of ~ 1-2.5 M0, that there is a dip in 
the trend. Whether this is due to the physical mass of the 
stars, or is an observational effect from different sample is 
unclear. This dip will be discussed further in Section [6.41 in 
regard to the HAeBes of this sample. An investigation into 
the meaning of this dip, and how the relationships behave 
in CTTs and HAeBes, is presented in a dedica ted paper by 
this group to the topic jMendieutfa et al .l2015h . Overall, the 
figure shows a trend that covers a large mass range spanning 
low mass CTTs to high mass HBes, with only some slight 
deviation in the HAe mass range. 

Figure [T3 instead shows a relationship of the luminosi¬ 
ties instead of the mass; specifically of how Lace changes 
as a function of L*. Again, comparisons are made against 
HAeBes and CTTs from the literature. A positive correla¬ 
tion between the two is also seen here of Lace oc 
This trend in the data shows a scatter of around 2 dex in 
Lace throughout the luminosity range covered; this scatter 
is comparable to the scatter in Macc shown in Figure [TT] 

In total, accretion rates, and therefore accretion lumi¬ 
nosities, have been calculated for 81 stars in the sample. 
Their values are seen to agree with previous literature esti¬ 
mates of accretion in HAeBes. The accretion rates obtained 
are observed to increase with both temperature and lumi¬ 
nosity; this trend is seen in the literature for CTTs and 
HAeBes alike. 

6 DISCUSSION 
6.1 Overall Results 


ADs is clearly detected in 62 of the stars, while a further 
26 stars have upper limits placed on them. The remaining 3 
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Figure 14. The left panel shows all the stars for which log(g) and Tgff could be determined from the spectra. These are translated into 
an HR-diagram in the right hand panel. In both panels the colour of the points reflect the strength of the accretion rate determined in 
each star; dark-blue symbols are the strongest accretors while light-red ones are the wea kest accretors. The squares d enote objects where 
Macc is an upper limit. The PMS evolutionary mass tracks of iBressan et all l|2012l ') and iBernasconi &: Maeden l|l996lj are also plotted as 
solid grey lines, and are labelled according to mass. Stars which were moved onto the ZAMS are not included in this plot. 


stars are measured as having a negative or zero ADb^ The 
possible reasons for this for each star are now discussed: 
V590 Mon is observed to have ADb = 0 within the errors, 
which is acceptable as it may not be accreting; PDS 281 has 
been listed previously as a possible evolved star (| Vieira et al.l 
l2003h . as such the parameters derived in this work may be 
incorrect based on our assumptions, and if it is evolved it 
is unlikely to be accreting; HD 94509 has very narrow and 
deep absorption lines in its spectrum which suggest it is a 
super-giant star with a low log(g), as supp orted by past ob¬ 
servations (IStephenson fc Sanduleaklll97ll l. while such low 
values are not covered by our adopted model atmospheres. 
Futher investigation into the accretion properties of these 
3 stars through emission lines will be presented in a future 
paper by the authors, though their questionable nature as 
PMS objects should be noted. 

There are 7 objects for which the measured ADb 
value cannot be reproduced though magnetospheric accre¬ 
tion shock-modelling, using the method we adopt. This is 
because the appropriate ADb vs. Macc curve calculated 
for each of the stars, based on its stellar parameters, can¬ 
not reach the observed ADb before a 100% filling factor is 
achieved (see Figure[9]for the points at which a 100% filling 
factor is seen for different temperatures). Within this subset, 
3 stars have a very large ADb of > 0.85 (PDS 133, R Mon, 
and DG Cir), 3 have temperatures exceeding 20000 K (HD 
141926, HD 53367, and HD 305298), while the final star lies 
in between these two scenarios having a strong ADb value 
and is mid-B spectral type (HD 85567). These stars are all 
HBes. 

Additionally, there are 12 stars whose measured ADb 
are modelled by filling factors of greater than 25% of the 


stellar surface. T his is allowed, but it is an unusual occur¬ 
rence under MA (IValenti. Basri fc Johnslll993l : iLong et al.l 
[2^. A filling factor greater than 1 is an unphysical value, 
as it implies that the accretion column covers more than the 
total surface area of the star. This suggests that the MA sce¬ 
nario adopted here needs to be revised, or discarded, for the 
stars with unphysical filling factors. Caution should be ex¬ 
ercised when considering the Macc values of stars with high 
filling factors. This amounts to 9% of ADb detections being 
non-reproducible though the adopted MA shock-modelling, 
with a further 15% having unusually high filling factors. All 
of this gives a possible indication that MA may not be ap¬ 
plicable in all HAeBes; particularly for stars with a large 
ADb, or which have high temperatures i.e. the HBes. The 
remaining 76% can be fitted successfully within the context 
of MA. 

6.2 HR-diagram 

Using the spectra, log(g) could be determined for the ma¬ 
jority of the sample in addition to Tgfj; for these stars their 
stellar parameters were determined using PMS tracks. Their 
placement on these tracks confirms the young nature of these 
stars and is shown in Figure im in the left-hand panel, while 
the right-hand panel shows the corresponding HR-diagram. 
The stars which required revised distances to be calculated 
for them (see Section 13.311 are not included in Figure [TT] 
as their placement is artificial compared to the other stars. 
A large proportion of the sample are clustered between 2- 
3 Mq, which is likely caused by a combination of two effects. 
The first being that lower mass sources are more n umerous, 
as described by the initial mass function, IMF, (ISalneterl 
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Table 4. Compares the number of HAeBes 
found in different mass bins with the theoreti¬ 
cal IMF distribution 


Mass Bin 

Theoretical 

This Work 

1-2 Mq 

99 

3lt 

2-4 Mq 

39 

36 

4-8 Mq 

15 

11 

8-16 Mq 

6 

11 

> 16 Mq 

4 

2 


t This sample is focused on HAeBes and does 
not represent the 1-2 Mq bin well, as HAeBes 
are generally more massive. 


Il955h . The range of masses determined in this work agrees 
fairly well with a typical Salpeter-IMF distribution, partic¬ 
ularly when considering the selection criteria (the criteria 
skews our sample towards high mass objects, as these are 
the ones of greater interest in this work). Table [T] shows the 
comparison of the m ass distribution in this wo rk versus the 
distribution given in IZinnecker fc Yorkd (l2007ll for a typical 
IMF function. 

The second aspect, which may be contributing to the 
clustering, is a visibility effect due to low mass stars be¬ 
ing more evolved and less extinct than younger high mass 
stars (as predicted by a comparison between the Kelvin- 
Helmholtz time-scale and the free-fall time-scale). This sec¬ 
ond point is supported by the Av values measured in this 
sample where, in general, the lower mass objects tend to 
have lower Av values. However, it should be noted that a 
high Av does not necessarily mean that the HAeBe has a 
high mass, as many low mass stars of young ages also have 
high extinction values (e.g. V599 Ori has an Av = 4.65, 
but only has Teff = 8000 K and M* = 2.5 Mq). An¬ 
other point to note is that clustering of the stars in the 
figure could be attributed to the stars actually belonging to 
the same cluster. The main star forming regions in which 
some of the HAeBes in this work appear to be associ¬ 
ated with are the Orion-OBl, Mon- OBl, CMa-Rl and Sco- 
OB 2 regions llde Zeeuw et al. 1 1 9 9'^: I Shevchenko et al.lll999l : 


iDahm fc Simonll2005l : Ivan LeeuwerJl^OTli . Since the regions 

are located at fixed distances, clustering of luminosities will 
occur if the stars are of similar spectral type. It is worth not¬ 
ing that the number of stars in each mass bin of a cluster is 
governed by the initial mass of the cloud in which they form. 
By looking at just a few star forming regions we naturally 
get clusters of similar mass stars in each one; resulting in 
clustered regions in an HR-diagram. However, only a small 
number of distances are adopted from the literature as an 
input parameter in this work, and they are drawn from var¬ 
ious catalogues and regions on the sky. The spread on the 
HR-diagram can simply be attributed to relatively low num¬ 
ber statistics. 


6.3 Age 

Next, the ages of the stars are investigated. Generally for 
PMS stars, the higher the mass, the younger the object is. 
Figure [TJ] shows age against temperature (Teff is roughly 
proportional to M* for MS stars, and stars close to the 
ZAMS). The plot shows an inverse relationship between age 



Figure 15. The PMS tracks in Figure [14] are used to obtain an 
age for each star. The ages determined are plotted here against 
the temperature of the star. The plot shows how the older stars 
are always the cooler stars i.e. the ones with a lower mass, which 
evolve towards the main sequence slower than their high mass 
counter parts, as expected. However, all of the hottest objects, 
T* > 20000 K, are seen to be the youngest ones, age< 0.5 Myr. 
Some cool and young stars are also present, which are likely to 
be in the early stages of their PMS evolution. 


and temperature where increasing temperature results in 
younger ages. This is as expected as the hot objects, which 
evolve faster, will move away from the PMS stage of their 
lives quicker than the lower mass stars, allowing the higher 
mass stars only to be seen at an early age. This point is 
worth stressing when it comes to looking at HAeBes statis¬ 
tically, as the HBes will always be much younger than the 
majority of HAes, but they can also be much closer to the 
main sequence, as this is relative to their mass. 

Figure [16] shows how Macc changes with the age of a 
star. As the age increases the accretion rate diminishes, 
much like what has been seen for the temperature. A fit to 
the data provides a relationship of Macc oc t”’’, were t is the 
age in Myrs, and r} = 1.92 ± 0.09. The hgure also shows the 
stars split into separate mass bins too. A fit to just the HBes, 
where M* > 3 Mq, obtains rj = 1.93 ± 0.24, which is very 
similar to the result for all of the HAeBes, only slightly off¬ 
set. However, for the HAes alone in this work a much steeper 
relationship is obtained of ?7 = 4.06 ± 0.53. The HAeBes 
as a whole, and the HBe case, agree with t h e HA eBes in¬ 
vestigated in the work of iMendiguti'a et akl (l2012l l , where 
the authors obtain 77 = 1.8. For C TTs a relationship ha s 
been observed where rj — 1.5-2 .8 llHartmann et al]|l998l l. 
This range also encompasses the case for the HAeBes as a 
whole, and the HBes. However, more recent studies suggest 
that the relatio nship for CTTs is actually lower than this , 
where rj = 1.2 ( ISicilia-Aguil ar. Henning fc Hartmannll2O10l : 
ICaratti o Garatti et al.ll^l 2 l l. which suggests that there is a 
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Figure 16. Plotted here are the derived accretion rates, from 
Table [3l against the age of the star, in a log-log plot. Various mass 
bins are noted by symbol and colour. Various fits are made to the 
different mass bins. All of the HAes are shown in red (< 3 Mq), 
and can be fitted by a relationshiip of Afacc cx t“4.06±0.53^ shown 
as a dotted red line; the remaining HBes (> 3 Mq) are all fitted 
with the dashed blue line where Macc ex t—1-9310.24 penally, a fit 
to all of the HAeBes is shown in black of Macc oe t“l-92±0.09 2 
CMa has been excluded from this fit (see the text for discussion). 

difTerence in the Macc oc relationship between the CTTs 
and the HAeBes. Some caution should be noted for the ages 
of the HBes as their ages are less accurate than the HAes, 
since the HBes are younger. In particular, there are some 
stars which are suspected to have ages < 1 Myr, the uncer¬ 
tainty in the ages of these stars are taken into account in 
the fitting. 

Overall, these relationships indicate that Macc could be 
an evolutionary property of HAeBes, which decreases as the 
star evolves; possibly accreting all of its material or dispers¬ 
ing its disc with time. Modelling of disc dispersion through 
photoevaporation suggests that the disc lifetimes are indee d 
shorter for more massive stars dGorti fc Hollenbachl 120091) . 
This offers an explanation for the steep exponent observed 
in the HAes, in which we could be observing the transition 
stage of disc dispersion as they approach the main sequence, 
resulting in a decreased accretion rate. The HBes, on the 
other hand, are younger and may not be dispersing their 
disc yet, which allows them to retain a more shallow rela¬ 
tionship between Macc and age. 

6.4 Accretion Rate vs. Stellar Parameters 

Moving on from the age, the next natural questions are: how 
is Macc related to the stellar parameters of the star; and are 
they influenced by it or vice-versa? In Figure [T7] a compar¬ 
ison is made between Macc and Tefj. The scatter in Macc 
remains constant at about 2 orders of a magnitude through¬ 



Figure 17. Macc is shown versus Tefj in a log-log plot here, 
where Macc appears to be increasing steadily with increasing 
temperature. The red squares denote upper limits. The outlier 
at log(Macc) ~ —3.0 and log(reff) ~ 3.95 is Z CMa. 

out. There is one object, Z CMa, which can be seen as an 
outlier from the general scatter. This star is cool, 8500 K, 
very massive, M* = 11 Mq, and has a very large Balmer 
Excess, ADs = 1.05. Its placement on the HR-diagram and 
PMS tracks puts it at a very early stage of evolution, in 
which it appears to be able to accrete at great rates. This 
star appears to be an exception to the majority of other 
stars and is excluded in all fitting because of this. The over¬ 
all trend is that Macc increases steadily with temperature; 
the temperature of a star is generally proportional to its 
mass leading to the next relationship. 

Figure [18] compares the log(Macc) vs. log(M*) rela¬ 
tionship, along with a series of fits to the data. In the 
hgure the stars are split into two groups, the HAes and 
the HBes, which comprise 60% and 40% respectively 
of the total sample (the split between the two regimes 
is made at 3 Mq). A best fit to the HAes is made of 
Macc oc while for the HBes a shallower re¬ 
lationship of Macc OC is seen. An overall ht 

to the HAeBes is obtained of Macc oc ^3.72±o.27^ which 
lies between the HAe and HBe regime and favours the 
HBe case, which covers a greater mass distribution. When 
considering the HAeBes as a whole, the relationship found 
between Macc and M* is a factor of ~2.0 larger in exponent 
than in low ma ss PM S stars, where Macc oc 
llMuzerolle et al.l l2005l : iNatta. Testi fc Randichl I 2 OO 6 I ) . 
The trend we observe of a steeper relationship between 
accretion r ate and stellar ma s s, over CTTs, agrees with the 
findings of iMendieutfa et al.l (l2011hl) . who also obtained a 
relationship with a much steeper exponent; they obtained 
Macc oc Although, the exponent for the HAeBes 

as a whole obtain e d in th is work is slightly shallower than 
iMendigutfa et al] (l201ld) . this could be attributed to 
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Figure 18. Macc versus M* is plotted here, for all stars in 
which Macc could be determined. The stars are split into HAes 
(as triangles, where M* < 3 Mq) and HBes (as squares where 
Mir > 3 Mq). Upper limits are denoted as the points in red with 
downward arrows from them. Separate fits are made to the HAes, 
HBes, and the group as a whole, of the form Macc oc M“, where a 
is found to be 8.59ibl.46, 2.82ib0.51, and 3.72ib0.27, respectively. 
A discussion of the fits is provided in the text. 


Figure 19. Lace versus L* is plotted here, for all stars in which 
Macc could be determined. The stars are split into HAes (as 
triangles, where M* < 3 Mq) and HBes (as squares where 
Mir > 3 Mq). Upper limits are denoted as the points in red 
with downward arrows from them. Separate fits are made to the 
HAes, HBes, and the group as a whole, of the form Lacc oc LJ, 
where a is found to be 1.15 zb 0.20, 0.85 zb 0.12, and 0.97 zb 0.06, 
respectively. A discussion of the fits is provided in the text. 


our sample containing more HBes. Overall, it is apparent 
that HAeBes have higher accretion rates and a steeper 
relationship to M* than CTTs. This could be due to 
HAeBes being younger stars, which are in earlier stages 
of accretion. The different relationships observed between 
HAes and HBes, could be due to the HAes crossing into 
a transitional disc phase, in which accretion rates may be 
lower (see the discussion above in Section lOl) . This is also 
possible due to the HAes having a longer stage of evolution 
compared to HBes. 

Alternatively, luminosities can be compared against 
each other instead of masses; Figure [T§] shows the lumi¬ 
nosity plot of log(I/acc/L 0 ) vs. log(I/*/L 0 ). A best fit to 
all the HAeBes is found of Lacc oc ij0.97±o.o6 'pjjjg jg 
in agreement with the work of iMendieutfa et al.l (l2011lj) . 
where these authors also found a shallower relationship 
for HAeBes of Lacc oc Ll'^. As in the previous compar¬ 
isons presented, when looking at the masses, the HAes and 
HBes are divided into two groups in the figure. A best 
fit to the HAes is obtained of Lacc oc 7 ;,i i 5 ±o. 2 o^ which is 
sligh tly shallower than the trends seen in CTTs of L acc oc 
Li ® (iNatta. Testi fc Randich|[2006l : lTilling et al.ll2008l i. The 
HBes, on the other hand, demonstrate an even shallower re¬ 
lationship of Lacc oc jg turn shifts the weight¬ 

ing when looking at the HAeBes as a whole. The data is 
suggestive of Lacc being tightly correlated with L*, but the 
exact relationship changes in exponent between the HAe and 
HBe regime. 


Table [S] is provided as a condensation of the various ac¬ 
cretion relationships extracted from the graphs. It demon¬ 
strates clearly the changes between the different mass groups 
when looking at accretion as either a function of mass, lu¬ 
minosity, or age. 


7 CONCLUSIONS AND FINAL REMARKS 

To conclude, we have presented the largest spectroscopic 
survey of HAeBes to date, obtaining the following results: 

• Basic stellar parameters are determined for the whole 
sample, in a homogeneous fashion, by initially constrain¬ 
ing Teff and log(g) using some of the best spectra available 
for these stars. As such, the parameters are more consistent 
between objects, while only a handful of objects require spe¬ 
cialist treatment. The findings are in agreement with previ¬ 
ous works in the literature. 

• A UV-excess, ADb, is clearly detected in 62 stars of 
the sample, with upper limits allowed on a further 26 stars. 
Macc is determined for 81 of these stars through modelling 
within the context of magnetospheric accretion. However, 7 
of the ADb detections cannot be reproduced in this context. 
These 7 stars are all HBe stars, often with very large ADb 
values of > 0.85, or high temperatures exceeding 20000 K. 
This suggests a possible breakdown in the MA regime for 
HBes, particularly for early-type HBes. 

• A clear trend is observed of Macc increasing as a func¬ 
tion of stellar mass. The relationship obtained is a power law 
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Table 5. Summarises the exponents describing accretion relationships for 
various PMS mass groups. The relationships are of the form A = , with 

C being the exponent. 


Group 

^acc 

VS. 

^acc 

-^'acc 

VS. 

La.cc 

Macc VS. age 

CTTs 

2.00 

± 

0.20 


1.50 

-1.20, -1.50 - -2.80 

CTTs-l-HAeBes 

2.47 

± 

0.07 

1.19 

± 

0.03 

N/A 

HAes 

8.59 

± 

1.46 

1.15 

± 

0.20 

-4.06 ± 0.53 

HBes 

2.82 

± 

0.51 

0.85 

± 

0.12 

-1.93 ± 0.24 

HAeBes 

3.72 

± 

0.27 

0.97 

± 

0.06 

-1.92 ± 0.09 


Refere nces for CTTs val ues _ are from: iMuzerolle et al] 

||200.4>: iNatta. Testi &: Randiclj ll2006h for Column 2 — 


iNatta. Testi Randichl 120061'); iTilline: et al.l i 

20081) for Column 3 - 

and 

Sicilia-Aeuilar. Hennine Sz, HartmannI ' 

2010 

): ICaratti o Garatti et al.l 

(12011 

) for the first value in Column 4, with 

Hartmann et al.l l|l998li for the 


second value. 


of the form: Macc oc ’^2=*=o.27 jg steeper law than 

previously observed in CTTs, which is only Macc oc . 
We interpret this increased exponent, in the relationship 
between Macc and M*, for HAeBes as a possible combina¬ 
tion of them being younger and therefore more active in 
formation than older CTTs. Deviations are seen between 
the HAes, where Macc oc jJBes, where 

Macc oc . An explanation could be that the HAes 

are crossing into a transitional disc phase, in which accretion 
rates may be lower. 

• There is also a trend between the accretion luminos¬ 

ity and the stellar luminosity, which is found to be Lace oc 
pio.97±o.06 sample. This is lower than found in 

CTTs where Lace oc L*® l|Natta. Testi fc R.andichl l200fil : 
iTilling et al.l[200^ ~l. However, for a subset of the HAes the 
relationship is much closer to the CTTs case, where we ob¬ 
serve Lace oc pi i5±o.20^ Contrast a shallower relation of 
Lace oc p0.85±o.i2 seen in the HBes. This demonstrates 
that the stellar luminosity of a star appears to be a good 
indicator of the accretion luminosity for a huge range of stel¬ 
lar luminosities, upto the HAe mass range, but there may 
be deviations in the HBe mass range. 

• A trend is also seen in the HAeBes between the age of 
the star and Macc, where the accretion rate decreases with 
increasing age, characterised by the form Macc oc with 
77 = 1.92 ± 0.09. This implies that the accretion rate de¬ 
creases as stars approach the main sequence. However, this 
result is affected by two factors. The first is the most massive 
stars, with the higher accretion rates ,are only observable at 
young ages due to their rapid evolution. Secondly, the less 
massive stars have a longer PMS lifetime, which could al¬ 
low their accretion rate to diminish within this time. These 
factors could explain the change in the relationship for the 
HAes case where r] = 4.06 ± 0.53. Overall, this suggests 
that the younger objects are indeed accreting at a faster 
rate, and that the accretion rate diminishes more quickly 
for older HAeBe stars, which could be a consequence of disc 
dissipation. 

This study has led to three main findings. Firstly, the 
HAeBes display relationships in accretion which are similar 
but different to CTTs. Macc is observed to have a steeper 
relationship with M* than seen in CTTs, while Lace shows 
a shallower relationship with L* than the CTTs case. Sec¬ 
ondly, there are also notable differences within the HAeBe 


group; when separating the HAes and HBes. Most notably, 
the HAes display a much steeper relationship in Macc when 
related to both age and M*. In both cases the steepness of 
the relationship is approximately double that seen in HBes. 
Although, the HAes also display a Lace relationship to L* 
which is comparable to CTTs. The third, and final, finding 
is that multiple early-type HBes, and stars with an observ¬ 
able ADfl of > 0.85, cannot be modelled successfully though 
magnetospheric accretion. This suggests that there is a pos¬ 
sible change in accretion mechanisms in these stars which 
requires further investigation. To further these findings, the 
next steps are to look at emission lines which are known 
tracers of accretion in CTTs and test their applicability in 
HAeBes. This will be presented in Paper H of the series. 
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APPENDIX A: PHOTOMETRY 

This appendix serves as a reference source for the sample. 
Provided here is Table El which contains all of the pho¬ 
tometry used in this work, along with references to the 
sources. Additionally, previously assigned literature values 
of distance and temperature are included for the whole 
sample (where possible, some do not have previous distance 
estimates). 


APPENDIX B: EXCEPTIONAL STARS 

Seven of the stars in the sample cannot be assigned a 
temperature from the spectra alone (see Section l3.ll) . For 
these objects a different approach must be undertaken on an 
individual basis in order to assign a limiting temperature. 
This is done by drawing upon as many literature sources 
on these objects as possible. Fortunately, there are very 
few objects in the sample which require this specialist 
treatment. The stars, and steps taken towards them are 
detailed below: 


VY Mon - This star is included here because it has 
the worst SNR of the sample. This makes accurate spectral 
typing difficult, but a cautious estimate of around 12000 K 
can be made for the tempera ture. This agrees with literature 
estimates of 8200-12000 K jMora et al.ri200ll : iManoi et al.l 
I 2 OO 6 II . A generous error of 4000 K is adopted. 

R Mon - In the spectra of R Mon all lines are seen 
in emission or as P-Cygni profiles, making any temperature 
estimate impossible from spectra alone. The temperature 
has been previously listed as around 12000 K in past works 
dMora et al.ll200i] : iManoi et al. lli^. We adopt this litera¬ 
ture temperature. 

Z CMa - Has lots of P-Gygni and emission lines in its 
spectra, but lacks absorption features for spectral typing. 
Again we must turn to the literature. In the literature this 
star is seen to have the large st spread in listed tempera¬ 
tures; ranging fro m 30000 K dvan den Ancker et al. 2004- 
Mmoi et al.ll200^, down to 11500 K dPonehew & BrittainI 


201lli and 8500 K dHinklev et al.l[2013l ') . We c hoose to adopt. 


and test, the most recent temperature from iHinklev et al.l 
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Table Al. Photometry from the literature. 


Name 

U 

(mag) 

B 

(mag) 

V 

(mag) 

R 

(mag) 

I 

(mag) 

Phot 

Ref 

Teff 

(K) 

Teff 

Ref 

D 

(pc) 

D 

ref 

UX Ori 

10.94 

10.71 

10.34 

10.12 

9.88 

a 

8410 

i 

340 

ap 

PDS 174 

13.54 

13.65 

12.84 

12.18 

11.42 

b 

18700 

b 

340 

ap 

V1012 Ori 

12.62 

12.46 

12.04 

11.61 

11.25 

c 

8600 

c 

340 

ap 

HD 34282 

10.15 

10.05 

9.89 

9.81 

9.71 

a 

8720 

i 

340 

ap 

HD 287823 

9.98 

9.90 

9.68 

9.59 

9.48 

b 

8720 

j 

340 

ap 

HD 287841 

10.63 

10.50 

10.21 

10.06 

9.89 

a 

8990 

i 

340 

ap 

HD 290409 

10.20 

10.11 

10.02 

9.96 

9.89 

b 

10500 

b 

340 

ap 

HD 35929 

8.71 

8.53 

8.12 

7.87 

7.61 

a 

6870 

k 

360 

aq 

HD 290500 

11.41 

11.35 

11.04 

- 

- 

d 

8970 

b 

470 

ap 

HD 244314 

10.42 

10.30 

10.10 

9.96 

9.80 

b 

8720 

1 

440 

ap 

HK Ori 

11.72 

11.79 

11.41 

11.05 

10.66 

a 

8460 

m 

440 

ap 

HD 244604 

9.68 

9.57 

9.38 

9.27 

9.12 

a 

8720 

1 

440 

ap 

UY Ori 

13.38 

13.16 

12.79 

12.56 

12.19 

b 

10500 

b 

510 

ap 

HD 245185 

10.02 

10.00 

9.91 

9.87 

9.82 

a 

9520 

1 

440 

ap 

T Ori 

11.38 

10.98 

10.43 

10.10 

9.63 

a 

8660 

i 

510 

ap 

V380 Ori 

10.80 

11.04 

10.53 

10.11 

9.50 

a 

9230 

n 

510 

ap 

HD 37258 

9.84 

9.80 

9.67 

9.59 

9.49 

a 

8970 

o 

510 

ap 

HD 290770 

9.18 

9.30 

9.27 

9.23 

9.18 

b 

10500 

b 

470 

ap 

BF Ori 

10.34 

10.05 

9.82 

9.68 

9.48 

a 

8990 

i 

510 

ap 

HD 37357 

9.00 

8.95 

8.84 

8.79 

8.72 

a 

9230 

1 

510 

ap 

HD 290764 

10.29 

10.20 

9.88 

9.68 

9.44 

b 

7200 

b 

470 

ap 

HD 37411 

10.07 

9.95 

9.82 

9.72 

9.58 

a 

9100 

1 

510 

ap 

V599 Ori 

17.07 

15.41 

13.76 

12.69 

11.55 

a 

7200 

b 

510 

ap 

V350 Ori 

11.39 

11.15 

10.82 

10.62 

10.34 

a 

8990 

i 

510 

ap 

HD 250550 

9.34 

9.61 

9.54 

9.32 

9.54 

e 

10750 

1 

280 

ar 

V791 Mon 

10.31 

10.68 

10.38 

10.12 

9.87 

b 

18700 

P 

1100 

as 

PDS 124 

13.15 

12.97 

12.44 

12.15 

11.81 

b 

9520 

b 

830 

w 

LkHa 339 

14.59 

14.24 

13.47 

12.80 

11.93 

a 

9230 

q 

830 

w 

VY Mon 

15.28 

14.56 

12.97 

11.82 

10.60 

a 

8200 

i 

800 

at 

R Mon 

12.17 

12.53 

11.93 

11.41 

10.87 

a 

12400 

i 

800 

at 

V590 Mon 

12.52 

12.75 

12.60 

12.42 

12.12 

a 

13000 

r 

800 

at 

PDS 24 

13.94 

13.62 

13.26 

12.98 

12.69 

b 

10500 

b 

590 

ap 

PDS 130 

14.42 

14.06 

13.40 

12.96 

12.44 

b 

10500 

b 

830 

w 

PDS 229N 

13.82 

13.70 

13.13 

12.74 

12.24 

b 

9520 

b 

830 

w 

GU CMa 

5.88 

6.56 

6.54 

6.47 

6.37 

a 

25000 

s 

1050 

au 

HT CMa 

12.55 

12.29 

11.87 

11.38 

11.87 

e 

9520 

q 

1050 

au 

Z CMa 

11.20 

10.50 

9.25 

8.40 

7.65 

f 

30000 

f 

1050 

au 

HU CMa 

11.72 

11.84 

11.55 

11.32 

11.16 

a 

11900 

q 

1050 

au 

HD 53367 

6.80 

7.37 

6.95 

6.67 

6.30 

a 

29500 

s 

1050 

au 

PDS 241 

12.33 

12.71 

12.06 

11.45 

11.11 

b 

30000 

b 

7000 

av 

NX Pup 

9.93 

9.96 

9.63 

9.38 

9.07 

a 

7290 

t 

410 

ap 

PDS 27 

14.61 

14.32 

13.00 

12.00 

10.98 

b 

17500 

u 

2900 

u 

PDS 133 

13.57 

13.61 

13.13 

12.79 

12.50 

b 

14000 

b 

2500 

b 

HD 59319 

7.86 

8.23 

8.31 

8.34 

8.42 

a 

11900 

V 

- 

- 

PDS 134 

12.50 

12.61 

12.20 

11.92 

11.65 

b 

14000 

b 

- 

- 

HD 68695 

10.00 

9.92 

9.82 

9.76 

9.66 

b 

9520 

w 

410 

ap 

HD 72106 

8.39 

8.50 

8.50 

8.49 

8.49 

b 

9810 

X 

370 

ap 

TYC 8581-2002-1 

12.18 

11.94 

11.48 

11.19 

10.95 

b 

t 8200 

b 

145 

ap 

PDS 33 

12.85 

12.63 

12.34 

12.16 

11.97 

b 

9520 

b 

370 

ap 

HD 76534 

7.68 

8.18 

8.07 

7.97 

7.84 

a 

20350 

y 

370 

ap 

PDS 281 

9.43 

9.46 

8.87 

8.50 

8.08 

b 

17050 

b 

370 

ap 

PDS 286 

14.39 

13.91 

12.15 

10.91 

9.76 

b 

30000 

b 

370 

ap 

PDS 297 

12.50 

12.34 

12.03 

11.83 

11.59 

b 

7850 

b 

145 

ap 

HD 85567 

8.11 

8.65 

8.51 

8.33 

8.08 

a 

12450 

z 

650 

aq 

HD 87403 

9.28 

9.31 

9.26 

9.22 

9.16 

b 

10100 

aa 

145 

ap 

PDS 37 

15.56 

15.06 

13.54 

12.38 

11.21 

b 

17500 

u 

3700 

u 

HD 305298 

10.36 

11.07 

10.86 

10.66 

10.47 

b 

36900 

ab 

- 

- 

HD 94509 

9.01 

9.15 

9.12 

9.10 

9.10 

a 

9730 

ac 

- 

- 

HD 95881 

8.53 

8.36 

8.19 

- 

- 

g 

8990 

ad 

118 

ap 

HD 96042 

7.89 

8.60 

8.47 

8.36 

8.23 

b 

25400 

ad 

- 

- 

HD 97048 

8.96 

8.80 

8.44 

8.20 

7.95 

a 

10010 

ae 

160 

aq 
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Table Al - continued 


Name 

U 

(mag) 

B 

(mag) 

V 

(mag) 

R 

(mag) 

I 

(mag) 

Phot 

Ref 

Teff 

(K) 

Teff 

Ref 

D 

(pc) 

D 

ref 

HD 98922 

6.74 

6.82 

6.77 

6.69 

6.61 

a 

10500 

w 

850 

aq 

HD 100453 

8.10 

8.07 

7.78 

7.60 

7.42 

b 

7390 

aa 

122 

aq 

HD 100546 

6.60 

6.70 

6.69 

6.67 

6.66 

a 

10500 

af 

97 

aq 

HD 101412 

9.57 

9.42 

9.24 

9.13 

9.00 

b 

10010 

aa 

118 

ap 

PDS 344 

13.09 

13.40 

13.15 

12.95 

12.77 

b 

15400 

b 

- 

- 

HD 104237 

6.64 

6.73 

6.52 

6.38 

6.23 

a 

8410 

z 

115 

aq 

V1028 Cen 

10.39 

10.70 

10.61 

10.48 

10.33 

a 

14100 

z 

130 

aq 

PDS 361S 

13.10 

13.35 

12.85 

12.49 

12.09 

b 

18700 

b 

- 

- 

HD 114981 

6.55 

7.13 

7.23 

7.27 

7.33 

b 

tl5400 

b 

550 

aq 

PDS 364 

13.85 

13.93 

13.46 

13.05 

12.63 

b 

11900 

ag 

118 

ap 

PDS 69 

9.92 

10.12 

9.80 

9.50 

9.12 

b 

17050 

ah 

630 

ah 

DG Cir 

15.96 

15.87 

14.75 

13.96 

13.06 

a 

15000 

ai 

700 

aw 

HD 132947 

8.87 

8.96 

8.91 

8.89 

8.89 

a 

10500 

af 

- 

- 

HD 135344B 

9.14 

9.14 

8.63 

8.16 

7.83 

h 

6590 

aj 

140 

ap 

HD 139614 

8.67 

8.64 

8.40 

8.26 

8.11 

b 

7850 

aj 

140 

ap 

PDS 144S 

13.59 

13.28 

12.79 

12.49 

12.16 

b 

8200 

b 

1000 

b 

HD 141569 

7.23 

7.20 

7.10 

7.03 

6.95 

a 

9520 

aj 

116 

aq 

HD 141926 

8.72 

9.20 

8.64 

8.21 

7.77 

b 

20300 

b 

- 

- 

HD 142666 

9.42 

9.17 

8.67 

8.35 

8.01 

b 

7580 

aj 

145 

ap 

HD 142527 

9.20 

9.15 

8.27 

- 

- 

g 

6260 

ak 

140 

ak 

HD 144432 

8.64 

8.53 

8.17 

7.94 

7.72 

a 

7350 

i 

160 

aq 

HD 144668 

7.28 

7.11 

6.78 

6.57 

6.38 

a 

7930 

al 

160 

aq 

HD 145718 

10.00 

9.62 

9.10 

8.79 

8.45 

b 

8200 

ag 

145 

ax 

PDS 415N 

13.43 

12.96 

12.04 

11.47 

10.85 

b 

7200 

b 

120 

ay 

HD 150193 

9.69 

9.33 

8.80 

8.41 

7.97 

a 

10010 

af 

120 

ay 

AK Sco 

9.56 

9.53 

8.90 

8.54 

8.18 

a 

6450 

am 

130 

aq 

PDS 431 

14.20 

13.99 

13.42 

13.02 

12.59 

b 

9520 

b 

145 

ap 

KK Oph 

13.17 

12.97 

12.36 

11.83 

11.03 

a 

8030 

an 

145 

ap 

HD 163296 

7.00 

6.96 

6.85 

6.80 

6.71 

a 

8720 

1 

119 

aq 

MWC 297 

14.94 

14.27 

12.03 

10.18 

8.80 

a 

23700 

ao 

250 

ao 


t These two stars are listed as objects QT3 (TYC 8581-2002-1) and QT4 (HD 114981) 
in the first Table of IVieira et alj ll2003ll . However, their places appear swapped in the sec¬ 
ond table by these authors. This swap is supported by additional photometry of HD 114981 
and by the authors observed temperatures and the temperatures derived in this work. Based 


on this we have sw apped the phot o metry fro m IV ieira et ^ 

stars. Refe rences: ( a 'l Ide Winter et alj l| 200l]). ('bl [Vieira et al.l _ _ _ 

lll999h. fd') iGuettei lll979lb fei iHerb st Shevchenkol lll99^. fP [v an den Ancker e t al.l ll2 004li. fg) 


ll2003l ') a round for these two 
2003h. I'd iMiros hnichen ko et al.l 


iMalf ait et alj ). fh') | C oulson j^Val^erl Jl9 9'^. fi') [Mora et alj ll200ll) , ^ (j) 

ll2005|i. fki IMiroshnichenko et al.l (1200^). Qi [Gray Corballyl l|l99'^ j 


Hernandez et al.l 

|Baine^^alJ j200^b 


n) I pinkenzeller ^^^to idtl lll98'^. f o) Icrav ^^ Corbahv lll993i . (ry) Cidale. Zorec &; Tringaniell3 
2 OO 1 I 1 . (q) iHer'nZnd^et al.l 1 I 2 OO 4 I'). (r) IPer ez^^aH l 20081 '). (s) iTiin A Dii^^al.1 ll200l | i. ('t ') 


fw') iHerbst &: Racind dlQTdh fx') iHouk 

1 I 1982 I). fv) fvalenti. Johns-Krull & Linskvl 1120001). (z) 

van den Ancker. de Winter & Tiin A Dih 

lll998l). (aa)lGuimaraes et alJ(l2006l). fab) Grahamll 

1970t). 

('ac') ISteohenson ^ SanduleakI lll97lh. ('ad'l 

iHouk Kr, GowlevI lll975l). fae) IWhittet et al. 

11987 

), (af) 
), (ai) 
), (al) 

Levenhaeen & Leisterl l2006h. fae) ICarmona et alj 1 I 2 OIOI'). (ah') iReiourth & Zinnecker 

1993 

Gahm ^ MalmortI lll98cl). fai^ iDunkin. Barlow RvanI iIIQQTI'). I'ak'l iFukaerawa et al, 

200e 

Tiin A Diie et al.1 lll989[). fam) lAndersen et al, 1 19891). fan) iHerbid ll2005l). fao) iDrew et alJ 1 

1997h. 



Dahm Sz Simo: 


(a aMvan I 
d ll2005h. 


(au) 


lll99Qh . fax') ^relbdsch &; Mamaiekl 


Shevchenko et ahl jlOOg^'Tfav^ jAvedisoval ||2000|~1 ~ 
20081 '). favl iLoinard et SI l|2008l b 


(aw) 


(I2013D . This is because their work spatially resolves the Her- 
big star in this system from its FU Or-like companion. In 
addition to this they provide SED fitting to the observed 
photometry to determine the temperature. 

PDS HI and PDS 37 - These two objects display very 
strong emission and P -Cygni profiles. They a re also the focus 
of a recent paper by lAbabakr et al.l ll2niRl . accepted) who 
determine distances and stellar parameters of the objects. 
We adopt their stellar parameters and distances in our work 


as they follow a similar methodology. The temperatures they 
found of ^21000 K, for both obj e cts, a re in agreement with 
the values found bv IVieira et al.l ll2003ll . 

PDS 133 - Another star devoid of any photospheric 
absorption in its spectra, and has extremely strong emission 
lines (the equivalent width of Ha is ~-100A). Therefore, 
using the spectra to assign a temperature is impossible. For 
this reason we ad opt a temperature around 14000 K, based 
on the literature llVieira et al.ll20()^ . 
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DG Cir - Another star with reasonably strong emission; 
the Balmer series are seen as P-Cygni profiles. A broad spec- 
tr al type of class B has previously been assigned to this sta r 
bv ISaniduleak fc Stephen^ and lVieira et ah l ll200^ . 

The authors Gahm fc MalmorT i igsd l do not give a spectral 
type but note its similarities to V380 Ori. A small indica¬ 
tion of absorption lines can be seen around 5200 A, but they 
appear close to many emission lines making an exact tem¬ 
perature determination difficult. We therefore agree with a 
B spectral type, and based on the absorption would narrow 
this to a late-B type star of -^11000 K, with a generous error 


of 3000 K. 













